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In a celebrated paper "Can one hear the shape of a drum?" M. Kac [Amer. Math. 
Monthly 73, 1 (1966)] asked his famous question about the existence of nonisometric 
billiards having the same spectrum of the Laplacian. This question was eventually 
answered positively in 1992 by the construction of noncongruent planar isospectral pairs. 
This review highlights mathematical and physical aspects of isospectrality. 
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I. INTRODUCTION 

Elastic plates are probably some of the oldest supports 
of sound production. They were used by most human 
cultures. Clay drums dated from the Chalcolithic have 
been found in graves in central Europe, and bronze drums 
dated from the second millenary B.C. have been dis- 
covered in Sweden and Hungary. However, it is usu- 
ally acknowledged that the scientific study of the vi- 
bration of elastic plates goes back only to the end of 
the 18th century, when the German researcher Ernst 
Chladni carried out the first systemati c investigations 



on the production of sound by plates ([Chladni 



1802 



Smilanskv and Stockmannl . 120071 ). When the plate was 
fixed in its middle and struck with a bow, it was set 
into vibration. The mode that was being excited was 
physically visualized by pouring sand on the plate: the 
sand accumulates at nodal lines, that is lines along which 
the plate does not oscillate. Some insight was brought 
into the mathematical theory of vibrating plates by the 
French mathematician Sophie Germain, who published 
Recherches sur la theorie des surfaces elastiques in 1821. 
In the course of the 19th century, Poisson, Kirchhoff, 
Lame, Mathieu, and Clebsch, devised analytic expres- 
sions for the description of the oscillation for elementary 
shapes such as the rectangle, the triangle, the circle, and 
the ellipse. 

The motivation for studying this problem was mainly 
that the wave phenomenon at the heart of membrane 
oscillations is in fact quite general. The stationary wave 
equation describing the problem arises in a variety of sit- 
uations. In many fields of physics, such as acoustics, 
seismology, hydrodynamics, and heat propagation, the 
mathematical formulation of the problem involves partial 
differential equations, and general solutions of these equa- 
tions can be found as superpositions of solutions of the 
so-called Helmholtz equation. In a d-dimensional space a 
stationary solution to the wave equation is an unknown 
function of d variables describing the problem, and the 
Helmholtz equation reads 



Af + Ef = 0, 



(1) 



where A is the d-dimensional Laplacian. Under suitable 
approximations, numerous problems can be cast in that 
form. For instance, in a certain regime the oscillations of 
the height / = f(x,y) of a thin vibrating plate at point 
(x,y) can be described by (J]). 

At the end of the 19th century, James C. Maxwell 
showed that the electric and the magnetic field behave 
like waves and established equations governing the time 
evolution of the electromagnetic field. From Maxwell's 
equations it is easy to prove that the electric and the 
magnetic field components also obey the same wave 
equation ([1} . Further interest developed in this equation 
when the wave-like behavior of matter was discovered in 
the early years of quantum mechanics. The Schrodinger 
equation was established in 1926 by Erwin Schrodinger 
to describe the spacetime evolution of a quantum sys- 
tem. The behavior of a particle can be described, in the 
framework of quantum mechanics, by a wave function ip, 
which is a function of the position of the particle, and 
which characterizes the probability amplitude il>(x) that 



the particle be located at a position x. If the system 
is described by the Hamiltonian H, the wave function 
satisfies the stationary Schrodinger equation = E^!, 
where E is the energy of the particle. For a particle 
of mass m and momentum p evolving in a box defined 
by its contour dB, the Hamiltonian describing the free 
motion inside the box reads H = p 2 /2m inside the box 
enclosure dB and oo outside, and the time-independent 
Schrodinger equation takes the form (JJJ). 

Mathematically, solutions of the Helmholtz equation are 
readily obtained in dimension d = 1. The problem of 
vibrating strings had been solved in the 18th century by 
Jean Le Rond d'Alembert. For a string of length L fixed 
at its two ends, solutions are simply given by f(x) — 
sm(mrx/ L), where n is an integer. The sound produced 
by the string has the possible frequencies nuo, with the 
fundamental frequency given by — c/(2L). 

Just as the one-dimensional case — which can describe a 
variety of physical situations — can be seen as a problem 
of vibrating strings, the two-dimensional case is usually 
studied from the perspective of its simplest mathematical 
equivalent, namely billiards. Billiards (in the mathemat- 
ical sense) are two-dimensional compact domains of the 
Euclidean plane R 2 . For instance, in quantum mechan- 
ics, the billiard models the behavior of a particle moving 
freely in a box whose dimensions are such that it can 
be approximated by a two-dimensional enclosure. The 
billiard problem is solved by looking for eigenfunctions 
ip and eigenvalues E that are solutions of Eq. {T]) inside 
the billiard, imposing boundary conditions on the bound- 
ary dB of the billiard. Physical problems impose specific 
boundary conditions. For instance hard wall domains in 
quantum mechanics impose that the wave function van- 
ishes on the boundary. In acoustics, clamping an elastic 
membrane imposes that the oscillations and their deriva- 
tive along the boundary vanish. The billiard problem 
usually considers the two following boundary conditions: 
Dirichlet boundary conditions ip\dB — 0, for which the 
function vanishes on the boundary, or Neumann bound- 
ary conditions 9 n V'|as — 0j f° r which the normal deriva- 
tive vanishes on the boundary. If such boundary condi- 
tions are imposed there is an infinite but countable num- 
ber of solutions to Eq. (fTJ) . We denote eigenfunctions of 
the operator —A by ip n and eigenvalues by E n , n e N, 
with < E\ < E% < E% ■ ■ ■ . Of course any combination 
of the above boundary conditions yields a different spec- 
tral problem. In this review however, we will be mainly 
concerned with Dirichlet boundary conditions. 

In the second half of the 20th century, quantum billiards 
were studied in the framework of quantum chaos. 
Quantum properties of classical systems were inves- 
tigated, and different behaviors were found according 
to the properties of integrability or chaoticity of the 
underlying classical dynamics. This quantum-classical 
correspo ndence led to various con jectures for integrable 
systems dBerry and Tabor . 1977 ) and chaotic systems 
( Bohigas et all , Il984f ). These conjectures rest on pow- 
erful mathematical tools that allow insight into the 
properties of solutions of the Helmholtz equation (JT]). 
For instance, the Weyl formula (see section IV. Aj) . or 
semiclassical trace formulas (see section IV.B.2[) , provide 
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a connection between the density of energy levels and 
classical features of the domains such as area, perimeter 
or properties of classical trajectories in the domain. 
The existence of such formulas and the conjectures on 
the quantum-classical correspondence indicate that the 
spectrum of a billiard contains a certain amount of 
information about the shape of the billiard. Therefore 
it is natural to ask how much information about the 
billiard can be retrieved from knowledge of the eigenvalue 
spectrum. For rectangular or triangular billiards, it 
is known that a finite number of eigenvalues suffices 
to entirely specify the sh ape of the billiard (see e.g. 
( Chang and Deturck . 19891 )). but is this true for more 
complicated shapes? 



In 1966, in a celebrated paper (|Kad . ll9 66), Mark Kac for- 
mulated the famous question "Can one hear the shape of 
a drum?" . This provocative question is of course to be un- 
derstood mathematically as follows: Is it possible to find 
two (or more) non-isometric Euclidean simply connected 
domains for which the sets {E n \\ n £ N} of solutions of 
dTJ) with * | Boundary = ^ are identical? More broadly, 
the question raises the issue of the inverse problem of 
retrieving information about a drum from knowledge of 
its spectral properties. As the spectroscopist A. Schuster 
put it in an 1882 report to the British Association for 
the Advancement of Science: "To find out the different 
tunes sent out by a vibrating system is a problem which 
may or may not be solvable in certain special cases, but 
it would baffle the most skillful mathematicians to solve 
the inverse problem and to find out the shape of a bell by 
means of the sounds which it is capable of sending out. 
And this is the problem which ultimately spectroscopy 
hopes to solve in the case of light. In the meantime 
we must welcome with d elight even the smallest step i n 
the desired direction." (jMehra and Rechenberg . l2000h . 
Actually, it was known very early, from Weyl's formula, 
that one can "hear" the area of a dru m and the length 
of its perimeter (see section IV. A| and ( Vaa et al . 2005h 
for a historical account of the problem). But could the 
shape itself be retrieved from the spectrum? That is, 
what kind of information on the geometry is it possible to 
gather from the knowledge of the spectrum, for instance, 
using semiclassical methods that allow investigation of 
the quantum-classical correspondence? And what kind 
of sufficient conditions allow the geometry to be entirely 
specified from the spectrum? 

Formally, an answer "no" to Kac's question amounts to 
finding isospectral billiards, that is non-isometric billiards 
having exactly the same eige nvalue spec trum. Since the 
appearance of Kac's paper ( Kac . 19661 ). far more than 
500 papers have been written on the subject, and innu- 
merable variations on "hearing the shape of something" 
can be found in the literature. Early examples of flat 
tori sharing the same eigenvalue spectrum were found 
in 1964 by Milnor in M 16 from nonisometric lattices 
of rank 16 in K 16 (see section IIII[) . Other examples 
of isospectral Riemannian manifol ds wer e cons tructed 
later, for example on lens spaces (IlkedaL [l980) or on 
surfa ces with constant negative curvature (| Vigneras . 
1980). In 1982, H. Urakawa produced t he first example s 
of isospectral domains in W\ n > 4 (|UrakawaL Il982h . 



(These examples are also described by (|Protterl . Il987l ).) 
More specifically, it is proved that there exist domains C 
and C in the unit sphere §™ _1 in K™, n > 4, which are 
Dirichlet and Neumann isospectral but not congruent in 
g?i-i This existence follows from the observation that 
there are finite reflection groups W and W that act on 
the same Euclidean space R™, n > 4, for which the sets 
of exponents coincide, and the intersections (C and C) 
of their cham bers with are not co ngruent in § n_1 . 

Then work of iBerard and Bessonl (|l980j ) is applied. 
In the late 1980s, various other papers appeared, giving 
necessary conditions that any family of billiards s harin g 
the same spectrum s hould satisfy (( Melrosel 1983 ). 
(jOsgood et aZ.l . ll988al ). (|Osgood et all . ll988bl )). and nec- 
essary conditions given as inequ alities on the eigenvalues 
were reviewed in (ProtterJ, [1987). 



But it was almost 30 years after Kac's paper that the 
first example of two-dimensional billiards having exactly 
the same spectrum was finally exhibited in 1992. The 
pair was found by C. Gordon, D. Webb and S. Wolpert 
in their paper "Isospectral plane domains and s urfaces 



via Riemannian orbifolds" ([Gordon et all Il992al ) . They 



gave a no as a final answer to Kac's question, and as a 
reply to Kac's paper, they publish ed a paper titled "On e 
cannot hear the shape of a drum" ( Gordon et all . Il992b ) . 
The most popularized example is shown in Fig.[T] Crucial 





FIG. 1 Paradigmatic pair of isospectral billiards with seven 
half-square shaped base tiles. The dotted lines are just for the 
eye. 

for finding the example was a theorem by Sunada (see 
section IVII.3j) asserting that when two subgroups are 
"almost conjugate" in a group that acts by isometrics 
on a Riemannian manifold, the quotient manifolds are 
isospectral. In fact, the other examples which were con- 
structed after 1992 all used Sunada's method. Later, the 
so-called transplantation technique was used, giving an 
easier way for detecting isospectrality of planar billiards. 
Still, essentially only 17 families of examples that say 
no to Kac's question were constructed in a 40 year period. 

Since the literature on isospectrality is large, and covers 
a broad spectrum of mathematical topics, we have cho- 
sen here to put the focus on isospectral billiards, that 
is, two-dimensional isospectral domains of the Euclidean 
plane, with Dirichlet boundary conditions. It is worth 
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noting that simple examples of isospectral domains can 
be constructed in the case of mixed Dirichlet-Neumann 
bo undary conditions. S uch constructions were proposed 
by (jLevitin et al\ , \200w (see section lVIILA[) . We now re- 
view some results on related topics, to which we will not 
return in this paper. 

First we mention several fun damental result s on isospec- 
trality that will be omitted. IZelditchl (1998) proved that 
isospectral simple analytic surfaces of revolution are iso- 
metric. That is, he considered the moduli space 1Z of 
metrics of revolution (S 2 ,<?) with the following proper- 
ties. Suppose that there is an effective action of S 1 by 
isometries of (S 2 ,<?). The two fixed points are N and S. 
Denote by (r, 9) geodesic polar coordinates centered at 
N, with 9 = being some fixed meridian jm from N to 
S. The metric g can then be written as g = dr 2 + a(r)d9 2 , 
where a : [0, L] h-> R+ is defined by o(r) = \S r (N)\/(27r), 
with |SV(A0| the length of the distance circle of radius 
r centered at N. The properties now are as follows: (i) 
g is real analytic, (ii) a has precisely one critical point 
Tq e]0,L[, with a"(r ) < 0, corresponding to an equa- 
torial geodesic 7^;, and (iii) the nonlinear Poincare map 
V lE for 7b is of twist type. 

Denote by 1Z* C R the subset of metri c s wit h simple 
length spectra in the sense of (jZelditchl . 119981 ). Then 



Zelditch p roved that Spe c : K* i-> R^. is 1-1. Further- 



more, in (IZelditchl . Il999h - see also (IZelditchl . l2000h . 
Zelditch showed that real plane domains fl that (1) are 
simply connected and real analytic, (2) are Z2 x Z2- 
symmetric (i.e., have the symmetry of an ellipse), and (3) 
have at least one axis that is a nondegenerate bouncing 
ball orbit, the length of which has multiplicity 1 in the 
length spectrum Lsp(p.), a re indeed determ ined by their 
spectrum. In recent work IZelditchl (|2004al ) pursued his 
goal of eventually solving the inverse spectral problem 
for general real analytic plane domains. We will return 
to this issue in more detail in section IV. J[ 

Concerning the known counterexamples in the plane, it 
should be remarked that the constructed domains are 
not convex (see e . g. Ap pendix |A")) . The objective of 
Gordon and Webbl ((l994l ) is to exhibit pairs of convex 
domains in the hyperbolic plane H 2 that are both Dirich- 
let and Neumann isospectral. They are obtained from 
nonconvex examples in the real plane by modifying the 
shape of a fundamental tile. Other interesting variations 
on the problem include the construction of a pair of 
isospectral (nonisometric) compact three-manifolds, 
called "Te tra" and "Didi", which h ave different closed 
geodesies ( Dovle and Rossettil 2004 ). 



The related question of graph isospectrality has also 
attracted much interest. We mention here a few results. 
A quantum graph is a metric graph equipped with a 
differential operator (typically the negative Laplacian) 
and homogeneous differential boundary conditions at 
the vertices. (Recall that a metric graph is a graph 
such that to each edge e is assigned a finite (strictly 
positive) length £ e £ R, so that it can be identified with 
the closed interval [0, £ e ] C R. Without the boundary 
conditions, the graph "consists of" edges with functions 
defined separately on each edge.) So there is a natural 



spectral theory associated with quantum graphs. Many 
results exist, and we just mention a few striking ones. 
One of t he main results in that spect ral theory can be 
found in ( Gutkin and Smilanskvl . l200ll ). where the trace 
formula is used to show that (under certain conditions) 
a quantum graph can be recovered from the spectrum of 
its Laplacian. (Necessary conditions include the graph 
being simple and the edges having rationally i ndependent 
lengths.) Using a spectral trace formula, iRothl ( 19841 ) 
in an e arly paper constr ucted isospectral quantum 
graphs. Ivon Below ( 2001 ). on the other hand, used 



the connection between spectra of discrete graphs and 
spectra of (equilateral) quantum graphs to transform 
isospectral discrete graphs i nto isospectral quantum 
graph s. Finally, we note that Parzanchevski and Band! 
(|2010r ) presented a method for constructing isospectral 
quantum graphs, based on linear representations of 
finite groups. Note that a di fferent notion of graph 
isospectrality was considered by iThasI (j2007bl ) based on 
the spectrum of the adjacency matrix of the graph. We 
return to this point in section [IV. A. 41 

To end this section, we give a short description of the 
contents of the paper. 

To familiarize the reader with the notions involved, we 
start by presenting a simple proof of isospe ctrality for the 
seminal example of Gordon et al. ( 1992al ) in section [TTJ 
Then the first historical examples of higher-dimensional 
isospectral pairs of flat tori are constructed (section ITlT|) . 
(Much more work has been done on isospectrality for 
the Laplace-Beltrami operator on flat tori in higher di- 
mensions than just the material we cover in section II. 
We refer to that section for more commentaries on that 
matter.) Section ITVl is devoted to the mathematical as- 
pects lying behind the construction of the known exam- 
ples of isospectral pairs. Then we review various aspects 
of the properties of isospectral pairs (section [Vj) , as well 
as experimental implementations and numerical checks 
of isospectrality (section IVIj) . As the first examples of 
isospectral billiards were produced by applying Sunada 
theory, a review of this theory is given in section IVIII 
In the last section we examine questions related to Kac's 
problem. 



II. A PEDESTRIAN PROOF OF ISOSPECTRALITY 

The first examples of isospectral billiards in the Eu- 
clidean plane were constructed using powerful mathemat- 
ical tools. We postpone these historical constructions to 
section lyiLEI The present section aims at illustrating the 
main ideas involved in isospectrality, so that the reader 
can acquire some intuition about it. More rigorous math- 
ematical grounds will be provided in the next sections. 



A. Paper-folding proof 

We start wit h a simple const ruction method that was 
proposed by IChapmanl (1995). It is based on the so- 
called " paper - folding" method. To illustrate it we follow 
(Thain, 2004j), where the method is illustrated on a simple 
example. 
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Consider the two billiards in Fig.[2j Each billiard is made 
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FIG. 2 The pair 73 (see Appendix |7\]) of isospectral billiards 
with a rectangular base shape. 

of seven identical rectangular building blocks. The solid 
lines are hard wall boundaries, the dotted lines are just a 
guide to the eye marking the building blocks. Let (j> be an 
eigenfunction of the left billiard with eigenvalue E. The 
goal is to construct an eigenfunction of the right billiard 
with the same eigenvalue, that is a function which: 

• verifies the Helmholtz equation (JXJ) ; 

• vanishes on the boundary of the billiard; 

• has a continuous normal derivative inside the bil- 
liard. 

The idea is to define a function ip on the right billiard as 
a superposition of translations of the function 4>. Since 
the Helmholtz equation ([1} satisfied by <p is linear, any 
linear combination of translations of (j> will be a solution 
of the Helmholtz equation with the same eigenvalue E in 
the interior of each building block of the second billiard. 
The problem reduces to finding a linear combination that 
vanishes on the boundary and has the correct continuity 
properties inside the billiard. The paper-folding method 
allows to satisfy all these conditions simultaneously. 



Take three copies of the left billiard of Fig. [5J Fold each 
copy in a different way, as shown in Fig. [3] (left column). 
Then the three-times folded billiards are stacked on top of 
each other as indicated in the right column of Fig. [3] note 
that the first shape (folding 1) has been translated on 
the left before being stacked, and that the second shape 
(folding 2) has been rotated by it in the plane of the 
figure. Once superposed, these three billiard yield the 
shape on the bottom right, which is the right billiard of 

Fig. m 

Now we make a correspondence between stacking two 
sheets of paper and adding the functions defined on these 
sheets; moreover, stacking the reverse of a sheet corre- 
sponds to assigning a minus sign to the function. For 
instance, in folding 3, a minus sign is associated in the 
right column with tiles 3 and 4, since they are folded back, 
and a plus sign is assigned to the other tiles since they 
are not folded. The function ip is defined by this "folding 
and stacking" procedure. For instance it is defined in the 
tile numbered 1 in the right billiard of Fig. [2] by 



■01 tile 1 = "4>\ti\c 1 + 4>\ti\c 2 



' tile 5- 



(2) 



The procedure above ensures that ip vanishes on the 
boundary and has a continuous derivative across the tile 
boundaries. Indeed, consider for instance the leftmost 
vertical boundary of the right billiard (i.e. the left 
edge of tile 1). On this boundary we have (p\ts\ e 5 = 
(since it is at the boundary of the left billiard), and 



0|tiio l = 4>\ti\e 2 since tiles 1 and 2 are glued together. 
Thus, ip given by Eq. ([2]) indeed vanishes on the leftmost 
vertical boundary of the right billiard. After we have 
checked by inspection all (inner and outer) boundaries, 
we have proved that the two billiards of Fig. [2] are 
isospectral. 



With the paper-folding method, it is clear that what mat- 
ters is the way the building blocks (the elementary rectan- 
gles in our example) are glued to each other, irrespective 
of their shape. We now show how the paper-folding proof 
generalizes to other shapes. Suppose we denote by 1, 2, 
and 3 respectively the left, right, and bottom edge of tile 
4 in the left billiard of Fig. [2j To obtain the whole billiard 
one unfolds tile 4 with respect to its side number 3, get- 
ting tile 7. Then tile 7 is unfolded with respect to its side 
number 2, yielding tile 6, and so on. The unfolding rules 
can be summed up in a graph specifying the way we un- 
fold the building block. The graphs in Fig@] correspond 
to the unfoldings yielding the billiards of Fig. [5] when ap- 
plied to a rectangular building block. The vertices of the 
graph represent the building blocks, and the edges of the 
graph are "colored" according to the unfolding rule, that 
is, depending on which of its sides the building block is 
unfolded. The graphs can alternatively be encoded by 
permutations a^\b^\ 1 < fi < 3. For instance for the 
first graph we have a W = (23) (56), = (12) (67), and 
a ( 3 ) = (25)(47). In fact, only three sides of the rectangle 
are involved in the unfolding. So we can start with any 
triangular-shaped building block, and unfold it with re- 
spect to its sides just as the billiards in Fig.[5]are obtained 
from the rectangular building block. This leads to billiard 
pairs whose isospectrality is granted by the paper-folding 
proof given above. For example, starting from the trian- 
gle in Fig. |4] and following the same unfolding rules, we 
get the pair of isospectral billiards shown in Fig. @] right. 
Taking a building block in the form of a half-square, we 
recover the example of Fig. [T] when the same unfolding 
rules are applied. 

The building block is in fact not even required to be a tri- 
angle or a rectangle. Any building block possessing three 
edges around which to unfold leads to a different pair of 
isospectral billiards. Another interesting example is ob- 
tained by taking a heptagon and unfolding it with respect 
to three of its sides, following the unfolding rules of Fig. [U 



This yi elds the first example produced by lGordon et al 
(|l992al lbh (see Fig.l5i 



Chapma produced more involved examples, fol 

lowing the same procedure. Starting from the build- 
ing block of Fig. [6] left, one obtains an ex ample of a 



pair o f chaotic billiards with holes. Similarly iDhar et al 



(|2003l ) constructed chaotic isospectral billiards based on 
the same idea: scattering circular disks were added in- 
side the base triangular shape in a way consistent with 
the unfolding. 

The central building block of Fig. [5] yields a simple dis- 
connected pair where each billiard consists of a disjoint 
rectangle and triangle. In this case, isospectrality can 
be checked directly by calculating the eigenvalues, since 
the eigenvalue problem can be solved exactly for triangles 
and half-squares. 




Slcc man and Hua (2000) considered a building block 
with piecewise fractal boundary: starting from a 
(tt/2, 7r/3, 7r/6) base triangle they cut each side into three 
pieces and remove the three triangular corners. Along the 
freshly made cuts a Koch curve is constructed, while the 





FIG. 4 Graphs corresponding to a pair of isospectral billiards: 
If we label the sides of the triangle by fi = 1 , 2, 3, the unfolding 
rule by symmetry with respect to side /j, can be represented 
by edges made of fi braids in the graph. From a given pair of 
graphs, one can construct infinitely many pairs of isospectral 
billiards by applying the unfolding rules to any shape. 



untouched sides still allow the Chapman unfolding (Fig. [5] 
right). This yields a pair of isospectral billiards with frac- 
tal boundary of dimension In 4/ In 3. 





FIG. 6 Examples of building blocks yielding isospectral pairs. 

A separate problem that will not be presented here 
is to find inhomogeneous vibrating membranes isospec- 
tral to a ho mogeneou s mem brane with the same shape 
(see, e.g., (iGottliebl. l2004h for circular membranes) . 



Knowles and McCarthy! (|2004l) used the isospectrality of 
the billiards of Fig. [T] to construct a pair of isospectral 
circular membranes by a conformal mapping. 



B. Transplantation proof 



FIG. 



5 Isospectral billiards. The top left figure is the seven- 



edged building block. From (|Gordon et all 1992a). 



The paper-folding proof can be made more formal 
be means of the so-called "tra nsplant a tion" method. 
This metho d was introduced in Berardl ( 1989 h_ Berard 
(Il99l Il993l). and discus sed by iBuser et all (| 19941) and 
Okada and Shudol (|200lh . It will be presented in more 
detail in section HVl Here we sketch the main ideas using 
a simple example. 

Consider the isospectral pair of Fig. [5J Let <f> be an eigen- 
state of the first billiard. Any point in the billiard can 
be specified by its coordinates a = {x, y) inside a build- 
ing block, and a number i arbitrarily associated with the 
building block (for example 1 < i < 7 in our example of 
Fig. [5]). Thus 4> is a function of the variable (a, i). Ac- 
cording to the paper-folding proof, a building block i of 
the second billiard is constructed from a superposition of 



three building blocks j obtained by folding the first bil- 
liard. We can code the result of the folding-and-stacking 
procedure in a matrix T, as 



T = 



( 


-1 


1 








1 





\ 




1 





1 








-1 
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-1 











-1 
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-1 


1 
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-1 


1 













1 





1 





1 
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1 








-1 








1 / 



(3) 



The paper-folding proof consists in showing that one can 
construct an eigenstate -0 of the second billiard as 



(4) 



where M is some normalization factor. That is, one can 
" transplant" the eigenfunction of the first billiard to the 
second one. The matrix T is called a "transplantation 
matrix" . The proof of isospectrality reduces to checking 
that ijj given by ©-(31) vanishes on the boundary and 
has a continuous derivative inside the billiard. 

Let us first transform the problem into an equiva- 
lent one on translation surfaces. Translation surfaces 
( Gutkin and Judge , 2000l) . also called planar structures, 
are manifolds of zero cu rvature w i th a finite number 
of singular points (see (jVorobetsl . Il996l ) for a more 
rigorous mathematica l defin ition). A construction by 
Zemlvakov and Katokl ( 1976h allows to construct a planar 
structure on rational polygonal billiards, that is polygo- 
nal billiards whose angles at the vertices are of the form 
on = nrrii/ni, with m^ni positive integers. This planar 
structure is obtained by "unfolding" the polygon, that is 
by gluing to the initial polygon its images obtained by 
mirror reflection with respect to each of its sides, and re- 
peating this process on the images. For polygons with an- 
gles on = irmi/rii, this process terminates and 2n copies 
of the initial polygon are required, where n is the gcd of 
the »j. Identifying parallel sides, one gets a planar struc- 
ture of genus in general greater than 1. This structure 
has singular points corresponding to vertices of the ini- 
tial polygon where the angle on = nmi/rii is such that 
rrij 7^ 1. The gen us of the translation surf ace thus ob- 



tained is given by (jRichens and Berrvl 119811 ) 



2 ^ 



1 



(5) 



A very simple example of a translation surface is the 
flat torus, obtained by identifying the opposite sides of 
a square. Such a translation surface corresponds to four 
copies of a square billiard glued together. 
The billiards of Fig. [2] possess one 27r-angle, two 37r/2- 
angles and eight 7r/2-angles each. The translation sur- 
faces associated to these billiards are obtained by gluing 
together 2n = 4 copies of the billiards, yielding planar 
surfaces of genus 4. They are shown in Fig. [7J Opposite 
sides are identified (e.g. in the first surface, the left edge 
of tile 1 is identified with the right edge of tile 5). Each 
surface has four singular points. The symbols o and • 
represent a 67r-angle, while the x and * symbols denote 
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FIG. 7 The pair 7s of isospectral billiards with a rectangular 
base shape unfolded to a translation surface (i.e. flat billiard 
with opposite sides identified). 



a 87r angle. An example of a straight line drawn on the 
first surface is shown on Fig. [7j The eigenvalue problem 
on these surfaces is equivalent to the problem on the bil- 
liards. It is however simpler to handle since the transla- 
tion surfaces have no boundary. Thus, only the continuity 
properties of the eigenfunctions have to be checked. 
Each translation surface is tiled by seven rectangles. 
Again, any point on the surface can be specified by its 
coordinates (a, i). Each tile on the translation surface 
has six neighboring tiles, attached at its left, upper left, 
upper right, right, lower right and lower left edge, and 
numbered from 1 to 6 respectively. For instance tile 1 is 
surrounded by: tile 5 on its left edge, tile 6 on its right 
edge, tile 3 on its upper left edge, tile 1 itself on its up- 
per right edge (because of the identification of opposite 
sides), tile 3 on its lower left edge and tile 1 on its lower 
right edge. The way the tiles are glued together can be 
specified by permutation matrices A^ v \ I < v < 6, such 
that — 1 if and only if the edge number v of i glues 
tile i to tile j. For instance for the first translation sur- 
face, the matrix specifying which tile is on the right of 
which is 

/ I 0\ 
1 

1 
A {2) = 1 (6) 

1 
1 

\ 1 / 

(tile 6 is on the right of tile 1, therefore A^ \ — 1, and so 
on). In a similar way, matrices B^ u \ 1 < v < 6, can be 
defined for the second translation surface. Now suppose 
there exists a matrix T such that 



Via 1< v< 6, A {v) T = TB {v) . 



(7) 
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Then for any given eigenstate <fi of the first translation 
surface we can construct an eigenstate V for the second 
translation surface, defined by Eq. In order to prove 
isospectrality we only have to check for continuity prop- 
erties at each edge. Suppose tiles I and j are neighbors. 
This means that there exists a v, 1 < v < 6, such that 
A\j =1. To prove the continuity of tp between tiles i 
and j, we have to show that the quantity 

C = #M)-V>(a,j) (8) 

is equal to zero for all a belonging to the edge between i 
and j. By definition of v we have A^) — 1 if and only if 
k = j. Therefore 



III. FURTHER EXAMPLES IN HIGHER DIMENSIONS 



«a,j) = ^4V(a,4 

k 

and C is given by 

C = V>(a,i)-^4^(a,fc). 

k 

Using Eq. (g]), we get 



(9) 



(10) 



C = Tik<K&, k) -Uj2 A ik T kk'H^ k'). (11 



k,k' 



The sum over k on the right-hand side yields a term 
(A^T)ik' . According to the commutation relation ([7]), 
it is equal to (TB^i)^, which gives 



(a,fc) 



0(a,fc') 



(12) 



Now the continuity of the function <j) ensures that all 
the terms between parentheses vanish. Thus C = 0, 
and continuity of ip is proved. Continuity of partial 
derivatives is proved in the same way. 

The proof rests entirely on the fact that we assumed the 
existence of a transplantation matrix T satisfying the 
commutation properties ([7]). It turns out that such a 
matrix exists. One can check that given the matrix 

/ 1 1 1 \ 

10 10 1 

1 1 1 

10 110 , (13) 

10 10 10 

1110 

V 1 1 1 / 

are satisfied for all v, 



T = 



the commutation relations ([7]) 
1 < v < 6. Thus the proof of isospectrality is completed. 
We return in section IIVI on this transplantation proof of 
isospectrality. 

A natural question is to know how one can find a suitable 
matrix T and permutation matrices A^ v \ B^ u ' verifying 
all commutation equations ([Jj). Historically these matri- 
ces were obtained by the construction of Sunada triples, 
as will be explained in section lVII.3l In fact, it turns out 
that the matrix T is just the incidence matrix of the graph 
associated with a certain finite projective space (the Fano 
plane in our example), as will be explained in detail in 
section IIVI 



Milnorl (|l964l) showed that from two nonisometric lattices 
of rank 16 in R 16 discovered by IWitti (|194ll ). one can 
construct a pair of flat tori that have the same spectrum 
of eigenvalues (all relevant terms are defined below). 
In this section, we describe a simple criterion for the 
construction of nonisometric flat tori with the same 
eigenvalues for the Laplace operator, from certain 
lattices (which was used by Milnor for the particu- 
lar case mentioned above), and then we construct, 
for each integer n > 17, a pair of lattices of rank 
n in R™ that match the criterion. Furthermore, we 
describe results of S. Wolpert and M. Kneser on the 
moduli space of flat tori. An interesting survey pa- 
per focused on the (elementary) const ruction theor y 
of isospectral manifolds has been given bv lBrooksl ( 19881 ). 



A. Lattices and flat tori 

A lattice (that is, a discrete additive subgroup) can be 
prescribed as ATI 1 with A a fixed matrix. For example, 
set 



.4 



1 
1 1 



(14) 



then the lattice Al? consists of the points of the form 



o(l,l) +6(0,1), a,bel. 



(15) 



An n- dimensional (flat) torus T is R n factored by a lattice 
L = AZ n with A e GL(n,R). The torus is thus deter- 
mined by identifying points that differ by an element of 
the lattice. 

If we return to the planar example above, the torus 
topologically is a donut — one may see this by cutting 
out the parallelogram determined by (1,1) and (0,1), 
and then gluing opposite sides together. 

With A,B e GL(n,R) are associated the lattices AZ n 
and Bl n . The tori R n /AZ n and R n /BZ n , B £ 
GL(n,R), are isometric if and only if All 1 and B1 n are 
isometric by left multiplication by an element of 0(n, R). 
The matrices A and B are associated with the same lat- 
tice if and only if they are equivalent by multiplication on 
the right by an element of GL(re, Z). So the tori W 1 j All 1 
and R™ I BIT 1 are isometric if and only if A and B are 
equivalent in 



O(n,R)\GL(n,R)/GL(n,Z)0 (16) 
Here, 0(n,R) is the orthogonal group in n dimensions. 



Let H, K be subgroups of the group G. Then the space of double 
cosets H \ G/K consists of the subsets ("double cosets") of the 
form HgK, with g £ G. (It is clear that G can be partitioned 
in these double cosets, and each such double coset itself can be 
partitioned in right cosets of H, and also in left cosets of if.) So 
in H\G j K , x ~ y if and only if there are h £ H and k £ K such 
that hxk = y. 
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The metric structure of R n projects to T, and 
volume(T) = |detA|; T carries a Laplace operator 



A 



(17) 



which is just the projection of the Laplacian of 
R n . The lengths of closed geodesies of T are given by 
|| a || for a arbitrary in L, || • || being the Euclidean norm. 

Let P be a symmetric matrix that defines a quadratic 
form on R™. The spectrum of P is defined to be the 
sequence (with multiplicities) of values 7 = N T PN for 
N £ Z n . The sequence of squares of lengths of closed 
geodesies of M. n /AZ n is the spectrum of A T A — Q; the 
sequence of eigenvalues of the Laplacian is the spectrum 
of ■iTT 2 (A~ 1 )(A^ 1 ) T = 47T 2 Q~ 1 . The Jacobi inversion for- 
mula yields for positive r, 



exp(-4n 2 TN T Q- 1 N) 



volume(T) 
(4ttt)™/ 2 



E 



exp(-MQM). 



(18) 



This equation therefore relates the eigenvalue spectrum 
of the torus to its length spectrum. We will see in sec- 
tion [VTb73] other examples of this connection between the 
spectrum of the Laplacian and the length spectrum. 



B. Construction of examples 



If L is a lattice of R n , L* denotes its dual lattice, which 
consists of all y £ M. n for which (x,y) £ Z for all x £ L; 
here, (•, •) is the usual scalar product on 1" x R n . Clearly, 
(L*)* = L, and two lattices L and L' are isometric if and 
only if L* and L'* are. 

Recall that two flat tori of the form R n /Lj, i £ {1,2}, 
are isometric if and only if the lattices Li and L2 are 
isometric. The following theorem gives a criterion for 
isospectrality of flat tori. 

Theorem III.l Let Li and L2 be two nonisometric lat- 
tices of rank n in R n , n > 2, and suppose that for each 
r > in R, the ball of radius r about the origin contains 
the same number of points of Li and L2. Then the flat 
tori R™/LJ and R^/L^ are nonisometric while having the 
same spectrum for the Laplace operator. 

Proof. Suppose x 7^ is an element of Li of length a. 
Then there is an a' < a such that the ball of radius 
a' centered at contains all elements of Li with length 
strictly smaller than a (since Li is discrete). For any a' < 
a" < a, the ball of radius a" centered at contains that 
same number of elements. This ball contains as many 
elements of L2 as of Lj., and since the ball centered at 
with radius a contains strictly more elements of Li, it 
follows easily that L2 also contains vectors of length a. 
Each element z £ Lj, i £ {1,2}, determines an eigen- 
function f(x) — e 2,K ( x ' z ) 1 for the Laplace operator on 
R"/L*, with corresponding eigenvalue A = (2ir) 2 (z,z), 
so the number of eigenvalues less than or equal to (2irr) 2 
is equal to the number of points of Lj contained in the 
ball centered at with radius r. 



We conclude that R n /L^ and R^/L?; have the same 
spectrum of eigenvalues, while not being isometric. ■ 



Milnor's Construction ay 
lattices in R 16 (Witt, 



By using the Wit t non- 



19411 ). iMilnorl (|1964T> 



isometric lattices m 
essentially used the aforementioned criterion to construct 
the first example of nonisometric isospectral flat tori. 



16 and L 2 6 



Starting from these two nonisomet ric latt i ces L 
of rank 16 in R 16 as described in Witt! (1941), one can 
in fact construct examples of isospectral flat tori in 
for all n, n > 16, as follows. The lattic es L} 6 



satisfy the condition of Theorem IIII.ll (iWittl . Il94l[ p. 



and Lj 6 



324). Now embed R 6 in R 7 in the canonical way. De- 
note the coordinate axes of the latter by X\ , X 2 , ■ ■ ■ , X17, 
such that (X 1 ,X 2 , . . . , X 16 ) = R 16 . Suppose £ ^ is a 
vector on the Xi7-axis which has length strictly smaller 
than any non-zero vector of Li (and L2). Define two 
new lattices Lj 7 (of rank 17) generated by Lj 6 and £, 
i = 1,2. Since Xn _L R 16 , it follows easily that for 
any r > 0, the ball centered at the origin with radius 
r contains the same number of elements of L} 7 as of Lg 7 . 
One observes that these lattices are nonisometric. Thus, 
by Theorem IIII.ll we obtain two nonisometric fiat tori 
= 1,2, which have the same spectrum of 



(17 IT 17 



eigenvalues for the Laplace operator. 

Inductively, we can now define, in a similar way, the 

nonisometric lattices L™ and of rank n, n > 17, 

satisfying the condition of Theorem IIII.ll and thus 

leading to nonisometric flat tori R™/L™*, i — 1,2, which 

have the same spectrum of eigenvalues for the Laplace 

operator. 



C. The four-parameter family of Conway and Sloane 

Let A be a positive-definite lattice. The theta function of 
A is: 



6a(t) = 



m=0 



N„ 



(19) 



where Im(r) > 0, and N m is the number of vectors x £ A 
of norm m. Oa can be thought of as a formal power 
series in the indeterminate q, although sometimes one 
takes q — e l7TT for further investigation, with r a complex 
variable. In that case, 0a (t) is a holomorphic function 
of r for Im(r) > 0. 



IConwav and Sloanel (1992) construct a four-parameter 
family of pairs of four-dimensional lattices that are 
isospectral (equivalently, that have the same theta se- 
ries (|19p ). In a similar way as before, such lattice pairs 
yield isospectral flat tori. The main construction of 
(jConwav and Sloand . 119921 ) is given by the next result. 



Theorem III. 2 (Conway and Sloane, 1992) Let 

Soo, eo, e\, e2 be orthogonal vectors satisfying 



-'OO '-'OO 



a/12, e -e = 6/12, e r ei = c/12, e 2 -e 2 = d/12, 
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where a,b,c,d > 0, and let [w, x, y, z] denote the vector 
weoo +xe + yei + ze 2 . Let u± = [±3, -1, -1, -l],u = 
[l,±3,l,-l],u± = [l,-l,±3,l],«a = [1,1,-1, ±3]. 
Then the lattices L + (a,b,c,d) spanned by i>+ , Vq , , v 2 
and L~(a, b, c, d) spanned by , Vq , v± , are isospec- 
tral. 



In this section we have seen that is essentially "easy" to 
construct (nonisometric) isospectral flat tori. The Milnor 
example was exhibited in 1964. But it has taken about 
30 years to find counterexamples to Kac's question in the 
real plane . . . 



Some small values of a, b, c, d give examples which 



were first found by ISchiemaim ( 1990f ). Substituting 
(a,b,c,d) = (7,13,19, 49), one obtains the pair of 
Earnest and Nipd (|l99lh . 



D. The eigenvalue spectrum as moduli for flat tori 

We now discuss some interesting results on the eigenvalue 
spectrum for flat tori. We already saw that there exist 
nonisometric isospectral flat tori. A natural question is 
now how such tori are distributed. 

The following theorem gives an insight into this question 
by considering the case of a continuous family of isospec- 
tral flat tori. 



Theorem III.3 (Wolpertl (|l978h ) Let T s be a contin- 
uous family of isospectral tori defined for s G [0, 1]. Then 
the tori T s , s £ [0, 1], are isometric. 

An interesting r esult by M. Kneser is the following (see 
( Wolpertl , 19781 ) for a proof). It states that, given an 
eigenvalue spectrum of some torus, only a finite number 
of nonisometric tori can be isospectral to it. 

Theorem III. 4 (M. Kneser) The total number of 
nonisometric tori with a given eigenvalue spectrum is fi- 
nite. 

The following result is rather technical. Its main message 
is that given two tori R n /AZ n and W n /BZ n with the 
same eigenvalue spectrum, then either these two tori are 
isometric, or the quadratic forms {A T A) and (B T B) lie 
on a certain subvariety in the space of positive definite 
quadratic forms. A more precise statement is as follows. 
Denote the space of positive definite symmetric n x n- 
matrices by p(n,R), and observe that the map 



A e GL(n,R) i-> A T A e p(n,R) 



(20) 



determines a bijection from 0(n,R) \ GL(n,M) to 
p(n,R). Then the following theorem holds. 



Theorem III.5 (|Wolpertl (|l978T )) There is a properly 
discontinuous group G n acting on p(n, R) containing the 
transformation group induced by the GL(n, Z) action 



S^A[Z], 



(21) 



where S € p(n,R) and Z € GL(n, Z). Given P, S € 
p(n,R) with the same spectrum, either g(P) = S for 
some g € G n , or P, S € Q3„, where the latter is a subva- 
riety of p(n,R). Moreover, 

(i) QJ„ = {Q e p(n,R) || spec(Q) = spec(R), R e 
p(n,R) with R ^ g(Q) for all g S G n }, and 

(ii) 2J„ is the intersection o/p(rt,R) and a countable 
union of sub spaces o/R m for some m. 



IV. TRANSPLANTATION 

The aim of this section is to describe the idea of trans- 
plantation in a more mathematical way than in scctionim 
This concep t was presumably first introduced by iBerard 
(|l99lll993h . There is in fact a deep connection between 
transplantation theory and the mathematical field of fi- 
nite geometries. First we review some elementary facts 
about finite geometries. Application of these tools to 
transplantation theory sheds light on the reasons for the 
existence of isospectrality. 



A. Tiling 

1. Graphs and billiards by tiling 



In this section, we follow lOkada and Shudd p001). 



Tiling. All known isospectral billiards can be obtained 
by unfolding polygonal-shaped tiles. As the unfolding is 
done along only three sides of the polygon we can essen- 
tially consider triangles. We call such examples isospec- 
tral Euclidean Tl-domains. The known ones are listed in 
Appendix El The way the tiles are unfolded can be speci- 
fied by three permutation dx d-matrices 1 < /i < 3 
and d £ N, associated with the three sides of the triangle 
and defined in the following way: M^f = 1 if tiles i and 

j are glued by their side fi; — 1 if the side (i of tile 

i is the boundary of the billiard, and otherwise. The 
number of tiles is, of course, d. Call the matrices 
"adjacency matrices" . 

One can sum up the action of the M^' in a graph with 
colored edges: each copy of the base tile is associated 
with a vertex, and vertices i and j, i ^ j, are joined 
by an edge of color /i if and only if — I. In the 

same way, in the second member of the pair, the tiles 
are unfolded according to permutation matrices N^', 
1 < fj, < 3. We call such a colored graph an involution 
graph for reasons to be explained later in this section. 
An example of such graphs is given in Fig. 01 If D is 
a Euclidean Tl-domain with base tile a triangle, and 
9Jt = {MW II /i € {1,2,3}} is the set of associated 
permutation matrices (or, equivalently, the associated 
coloring), denote by r(D,9JT) the corresponding involu- 
tion graph. 



The fo llowing proposition is easy but rather useful dThasl . 
2007bl ). 



Proposition IV. 1 Let D be a Euclidean TL-domain with 
base tile a triangle, and let ffl = {M M \\ \i e {1,2,3}} 
be the set of associated permutation matrices. Then the 
matrix 



3 

A v = E 



M 



00 



00 



Stj) 



(22) 
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FIG. 8 Two isospectral billiards with a triangular base shape 
on seven tiles. 



where Sij is the Kronecker symbol, is the adjacency 
matrix of T(D,Tl). u 



Setting. Let / be an eigenfunction of the Lapla- 
cian with eigenvalue A € M for the Dirichlet problem 
corresponding to the left-hand billiard in Fig. [5] Let 
/o, fx, . . . , /g denote the functions obtained by restriction 
of / to each of the seven tiles of the left-hand billiard, 
as indicated on the left in Fig. [8j For the sake of 
convenience, we write i for /j. 

The Dirichlet boundary condition is that / must vanish 
on each boundary segment. This is equivalent to the 
assertion that / goes into —/if continued as a smooth 
eigenfunction across any boundary segment; in fact, it 
goes into / o a where a is the reflection on the boundary 
segment. 

On the right in Fig. [8j we show how to obtain from / 
another eigenfunction of eigenvalue A for the right-hand 
domain. We define the function 1 + 2+4 which is actually 
the function 



TVansplantability. Two billiards are said to be trans- 
plantable if there exists an invertible matrix T — the 
transplantation matrix — such that 



f\ °Tl +/ 2 °T 2 +/4O r4, 



(24) 



TM (p) = N M T for all 



(23) 



If the matrix T is a permutation matrix, the two domains 
would just have the same shape. One can show that 
transplantability implies isospectrality, as seen in section 



We n ow discuss an exa mple exhibited by | Buser et 
(|l994l ). and first found bv lGordon et all (jl992al) . 



2. The example of Gordon et al. 



Buser ( 19881) constructed a pairs of isospectral flat 
surfaces 9Jti and OT2 as covers of a certain surface 
9Jto, using a pair of a l most c onjugate subgroups of 
SL(3,2). iGordon et a/1 (|l992af ) similarly constructed 
orbifolds 0\ and O2, respectively being the quotient 
by an involutive isometry of 9Jli and | Ox and 2 
have a common orbifold cover — it is the quotient by an 
involutive isometry of the common cover of 9Jti and SOT2 • 
The Neumann orbifold spectrum of Oi is precisely the 
Neumann spectrum of the underlying manifold SDi(Oi), 
and these latter underlying spaces are simply connected 
real plane domains. Furthermore, Dirichlet isospectrality 
of 97l(Oi) and 9Jl(02) is obtained by exploiting the 
Dirichlet isospectrality of SCJti and 9H2 ■ 

We now analyze this pair of isospectral but non- 
congruent Euclidean do mains. We follow the very 
transparent approach of Buser et~al\ |l994) to show 



isospectrality. As the reader will notice, this will in fact 
be an easy approach to (and example of) transplantabil- 
ity. 



where for k = 1,2,4, is the isometry from the central 
triangle of the right-hand billiard to the triangle labeled k 
on the left-hand one. Now we see from the left-hand side 
that the functions 1,2,4 continue smoothly across dotted 
lines into copies of the functions 0, 5, —4 respectively, so 
that their sum continues into + 5 — 4 as shown. Simi- 
larly way one observes that this continues across a solid 
line to 4 — 5 — (its negative), and across a dashed line 
to 2 — 5 — 3, which continues across either a solid or dot- 
ted line to its own negative. These assertions, together 
with the similar ones obtained by cyclic permutation of 
the arms of the billiards, suffice to show that the trans- 
planted function is an eigenfunction of the eigenvalue A 
that vanishes along each boundary segment of the right- 
hand domain. 

We have defined a linear map which for each A transforms 
the A-eigenspace for the left-hand billiard into the A- 
eigenspace for the right-hand one. This is a non-singular 
map (the corresponding matrix is non-singular), and so 
the dimension of the eigenspace on the right-hand side 
is larger than or equal to the dimension on the left-hand 
side. By symmetry, it follows that the dimensions are 
equal. Since A was arbitrary, the two billiards are Dirich- 
let isospectral. 



3. The other known examples 

A similar techniq ue as in the previous section allowed 
iBuser et al. ( 19941 ) to show that the series of billiard pairs 
they produced are indeed isospectral. All these pairs are 
listed in Appendix [XJ they were first found by search- 
ing for suitable Sunada triples, and then verified to be 
isospectra l (in the plane) by the tr ansplantation method 
(see also (|Okada and Shudol . l200lh for a further discus- 
sion about the subject of this section). 



2 Orbifolds are generalizations of manifolds; they are locally mod- 
eled on quot ients of ope n subsets of M n by finite group actions. 
We refer to ( Scot]], Il983f) for a formal introduction. 



4. Euclidean Tl-domains and their involution graphs 

To conclude this section, we address a related problem, 
namely isospectrality of the involution graphs associated 
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with the isospectral billiards. We say that two (undi- 
rected) graphs are isospectral if their adjacency matrices 
have the same multiset of eigenvalues. Note that this def- 
inition of graph isospectr ality is different from the def i- 
nition introduced in e.g. (jGutkin and Smilanskvl . 120011 ). 
where the spectrum of a metric graph is defined as the 
spectrum of the Laplacian on the graph whose edges are 
assigned a given length. 

The following question was posed bv iThasI (|2007bl ): Let 
(Di^Di) be a pair of nonisometric isospectral Euclidean 
Tl-domains, and let r(£>i) = T{D 1 ,{M^ || fi <E 
{1,2,3}}) and T{D 2 ) = T(D 2 ,{N^ || » € {1,2,3}}) 
be the corresponding involution graphs. Are T(Di) and 
r(Z?2) isospectral? Note that one does not require the 
domains to be transplantable. (The term "cospectral- 
ity" is also sometimes used in graph theory, instead of 
"isospectrality".) 

We now show that the answer is "yes" when the do- 



mains are transplantable. The proof is taken from (jThas 
2007bh . 



Theorem IV. 2 Let (_Di,£>2) be a pair of nonisomet- 
ric isospectral Euclidean Tl-domains, and let T(D\) = 
r(fli,{MW || // G {1,2,3}}) and T(D 2 ) = 
T(D2,{N^ || /! 6 {1,2,3}}) be the corresponding in- 
volution graphs. Then T(Di) and r(£>2) are isospectral. 

Proof. Define, for /i = 1,2, 3, M* as the matrix which 
has the same entries as except on the diagonal, 

where it has only zeros. Define matrices analo- 
gously. Suppose that TM^T' 1 = N^) for all p. Note 
the following properties: 



r* M and Nl"\ h = 
matrices, with at most one 1 entry on each row; 



• M?> and N { /> , /i = 1,2,3, are symmetric (0,1)- 



(m)i 



r<X> 



[M, 

and [A& ] ] m = 1$, where 



4f] 



1 if there is a 



1 on the i-th row of M* , and otherwise, if m is 
even, // = 1, 2, 3, and similar properties hold for the 

• TriM^M^) = Tt(M^M^) = for i ^ j and 

Tr (N^N^) = Tr(Ni j) N^) = for i ^ j; 

• Tr{Mi i] m¥ ] Ml k) ) and Tr(iV* (i) JV» (fc) ) are inde- 
pendent of the permutation (ijk) of (123) (this is 
because the individual matrices are symmetric); 

• the value of all traces in the previous property is 
(note that, if {i,j, k} = {1,2,3}, such a trace 
equals since the existence of a nonzero diagonal 
entry of M^M^M^, respectively N^N^N^, 
implies T(Di), respectively r(Z3 2 ), to have closed 
circuits of length 3). 



From Proposition HVT] h follows that A = J2l=i is 
the adjacency matrix of T(Di), and B = Y^L=i N* , the 
adjacency matrix of T{D2). 

Consider a natural number neNo. Then, with the pre- 
vious properties in mind, it follows that 



Thus by^ the follo wing lemma (cf. 

( van Dam and Haemerst 20031 Lemma 1)) the adja- 
cency matrices of T(Di) and r(Z?2) have the same 
spectrum. 

Lemma IV. 3 Two kx k-matrices K and K' are isospec- 
tral if and only if Tr(iT') = Tr(K' 1 ) for I = 1, 2, . . . , k. 



In section IVIIII we will see that other graph theoretical 
problems turn up in Kac theory. 



B. Some projective geometry 

There is a fascinating relation between the structure of 
isospectral billiards and the geometry of vector spaces 
over finite fields. In section [TT| we constructed pairs of 
isospectral billiards using unfolding rules. These unfold- 
ing rules can be encoded into graphs, like the ones in 
Fig. |U Thus the structure of a pair of isospectral bil- 
liards is entirely encoded into a pair of graphs that have 
certain specific properties. The graphs of Fig. @] are "col- 
ored" according to a certain set of permutations. It turns 
out that the group generated by these permutations is 
precisely the automorphism group of a projective space 
over a finite field, the so-called Fano plane. The Fano 
plane has many beautiful properties and appears in vari- 
ous places, such as combinatorial problems or the multi- 
plication table of the octonions. A representation of this 
finite projective plane is given in Fig. [9] Here we will 
see that the adjacency matrix of the graph representing 
the Fano plane is nothing but the transplantation matrix 
between the two isospectral billiards of Fig. 0] 
In order to understand this deep connection, basic no- 
tions of finite geometries and design theory are required. 
In this section we provide the necessary tools. More 
det ails about t he no tions considered here can be found 
in dHirschfeldl . ll99Sh . Note that some remarks about 
iso spectrality, projective geom etry and groups are made 
in ( Vorobets and StepirJ . Il998h . however the results there 
are not fully mathematically rigorous. 



Tr(A n ) = Tr{B T 



(25) 



1. Finite projective geometry 

Let F g be the finite field with q elements, q a prime power, 
and denote by V(n, q) the n-dimensional vector space 
over F 9 , n a nonzero natural number. Define the (n — 1)- 
dimensional projective geometry PG(n — 1, q) over ¥ q as 
the set of all subspaces of V(n, q). Note that PG(n-l, q) 
is often called the "Desarguesian" or "classical" projec- 
tive space. The projective space PG(— l,q) is the empty 
set, and has dimension — 1. 

Points in PG(n, g) correspond to one-dimensional 
subspaces of V(n, q), lines in PG(n, q) correspond 
to two-dimensional subspaces of V(n,q), and so on. 
Any (i-dimensional subspace of PG(n, q) contains 
(q d+1 — l)/(q— 1) points. In particular, PG(n, q) itself 
has (<7 n+1 -l)/(<7-l) points. It also has (g" +1 -l)/(<?-l) 
hyperplanes (i.e. (n — l)-dimensional subspaces). 

Example. The Fano plane PG(2, 2) shown in Fig. M 
has seven points and seven hyperplanes or lines (one of 
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FIG. 9 The Fano plane. 

which is represented as a circle in Fig. [SJ . Any line con- 
tains three points (we say that three points are "incident" 
with each line) and any point belongs to three lines (we 
say that three lines are "incident" with each point). The 
use of the word "incident" in both cases enhances the 
symmetry between points and lines in this geometry. It 
is precisely this geometry that lies at the root of isospec- 
trality. 

2. Automorphism groups 

The automorphism groups of finite projective spaces play 
a key role in isospectrality as the generators of these 
groups allow us to construct the graphs that encode the 
unfolding rules for the billiard construction. We now de- 
fine these groups and mention some of their properties. 
For group theoretical notions not explained here, we refer 
to the beginning of section IVIII 

An automorphism or collineation of a finite projective 
space is a bijection of the points that preserves the type 
of each subspace (i.e. lines are mapped to lines, and more 
generally d-dimensional spaces to d-dimensional spaces) 
and preserves incidence properties (i.e. intersecting lines 
are transformed into intersecting lines, etc...). It can be 
shown that any automorphism of a PG(n, q), n > 3, 
necessarily has the following form: 



(26) 



where A E GL(n + 1, q), a is a field automorphism of ¥ q , 
the homogeneous coordinate x = (xq, x\, . . . , x n ) repre- 
sents a point of the space (which is determined up to a 
scalar), and x CT = (xq,x°, . . . ,x°) (recall that x° is the 
image of Xi under a) . 

The set of automorphisms of a projective space natu- 
rally forms a group, and in case of PG(n, q), n > 3, 
this group is denoted by PrL(n + l,q). The normal 
subgroup of PTL(n +l,q) which consists of all automor- 
phisms for which the companion field automorphism a 
is the identity, is the projective general linear group, and 
denoted by PGL(n + l,q). So PGL(n + 1, q) = GL{n + 
l,q)/Z(GL(n+l, q)), where Z(GL(n+l, q)) is the central 
subgroup of all scalar matrices of GL(n+ 1, q). Similarly 
one defines PSL(ra+l, £?) = SL(n+l, q)/Z(SL(n + l, q)), 
where Z(Sh(n+ 1, q)) is the central subgroup of all scalar 
matrices of SL(n + X,q) with unit determinant. Recall 
that SL(n + 1, q) consists of the elements of GL(n + 1, q) 
with unit determinant. 



An elation of PG(n, q) is an automorphism of which the 
fixed points structure precisely is a hyperplane, or the 
space itself. A homology either is the identity, or it is an 
automorphism that fixes a hyperplane pointwise, and one 
further point not contained in that hyperplane. 



3. Involutions in finite projective space 

Let PG(n,q), n E NU { — 1}, be the n-dimensional pro- 
jective space over the finite field ¥ q with q elements, so 



that q is a prime power; we have |PG(n, g)| — /- 1 ' 
(Note again that PG(— l,q) is the empty space.) 
We discuss the different types of involutions that can oc- 
cur in the a utomorphism group of a finite projective space 
(|SegreLll96lj> . 



• Baer Involutions. A Baer involution is an in- 
volution that is not contained in the linear auto- 
morphism group of the space so that q is a square, 
and it fixes an n-dimensional subspace over F ^ 
pointwise. 

• Linear Involutions in Even Characteristic. 
If q is even, and 9 is an involution that is not of 
Baer type, 9 must fix an m-dimensional subspace 
of PG(n,g) pointwise, with 1 < m < n < 2m + 1. 
In fact, to avoid trivialities, one assumes that m < 
n — 1. 

• Linear Involutions in Odd Characteristic. 
If 9 is a linear involution of PG(n, q), q odd, the set 
of fixed points is the union of two disjoint comple- 
mentary subspaces. Denote these by PG(k,q) and 
PG(n -k-l,q),k>n-k-l> -10 

We are now ready to explore a connection between Inci- 
dence Geometry and Kac Theory. 



C. Projective isospectral data 

1. Transplantation matrices, projective spaces and isospectral 
data 

Suppose one wants to construct a pair of isospectral bil- 
liards, starting fro m a planar polyg onal base shape. The 
idea described in ( Giraudl . l2005al) is to start from the 
transplantation matrix T, and choose it in such a way 
that the existence of commutation relations 



(27) 



for some permutation matrices AfW, N**' will be known 
a priori. This is the case if T is taken to be the incidence 
matrix of a finite projective space; the matrices M^' 
and are then permutations of the points and the 

hyperplanes of the finite projective space. 

An (AT, k, A)— symmetric balanced incomplete block design 
(SBIBD) is a rank 2 incidence geometry, defined on a set 



3 We do not consider the possibility of involutions without fixed 
points, as they are not relevant for our purpose. 



14 



of N points, each belonging to N subsets (called blocks) 
such that each block is incident with k points, any two 
distinct points are contained in exactly A blocks, and each 
point is incident with k different blocks. 
Example. The points and hyperplanes of an n- 
dimensional projective space PG(n,q) defined over ¥ q 
is an (N, k, A)-SBIBD with N = (q n+1 - \)/(q - 1), 
k = (q n - l)/(q - 1) and A - (q^ 1 - l)/(q - 1). 
So the Fano plane is a (7, 3, 1)-SB1BD. 

Let T be an (TV, k, A)-SBIBD. The points and the blocks 
can be labeled from to TV — 1. One can define an N x N- 
incidence matrix T describing to which block each point 
belongs. The entries Ty of the matrix are Tij = 1 if the 
point j belongs to the block i, and otherwise. It is easy 
to see that the matrix T verifies the relation 

TT T = AJ + (iV — k)X/(k — 1)1, (28) 

where J is the N x TV-matrix with all entries equal to 
1 and I the N x N identity matrix. In particular, the 
incidence matrix of PG(n,g) verifies 

TT T = AJJ + (k - A)I (29) 

with k and A as given above. 

Any permutation a of the points of a finite projective 
space can be written as a d x d permutation matrix M 
defined by M ia ^ = 1 and the other entries equal to zero. 
Here d is the number of points. If M is a permutation 
matrix associated with an automorphism of the space, 
then there exists a permutation matrix TV such that 

TM = NT. (30) 

In other words, (|30p means that permuting the columns 
of T (which correspond to the hyperplanes of the space) 
under M is in some sense equivalent to permuting the 
rows of T (corresponding to the points of the space) 
under N . The reason that this occurs is the concept 
of "duality" ; in a finite projective space the points and 
hyperplanes play the same role. 

Consider a finite projective space 7r = PG(n, q) with in- 
cidence matrix T. With each hyperplane in 7r we asso- 
ciate a tile in the first billiard, and to each point in 7r we 
associate a tile in the second billiard. If we choose per- 
mutations ikfw in PrL(n + 1, q), then the commutation 
relation (|30[) will ensure that there exist permutations 
JV'Cw verifying 

TM M = N^ l) T. (31) 

Since these commutation relations imply transplantabil- 
ity, they also imply isospectrality of the billiards 
constructed from the graphs corresponding to and 

Constraints. If the base tile has r sides, we need to 
choose r elements 1 < /j, < r, in PTL(n + l,q) in 

such a way that (at least) the following remarks are taken 
into account. 

• Since represents the reflection of a tile with 

respect to one of its sides, it has to be an involution. 



• In order that the billiards be connected, no point 
should be left out by the matrices M^' — in other 
words, the graph associated to the matrices M^) 
should be connected. 

• If we want the base tile to be of "any" shape, there 
should be no closed circuit in the graph (in other 
words, it should be a finite tree). 

Assume one is looking for a pair of isospectral billiards 
with d = (q 3 — l)/(q — 1) copies of a base tile having the 
shape of an r-gon, r > 3. We need to search for r invo- 
lutions such that the associated graph is connected and 
does not admit a closed circuit. Such a graph connects 
d vertices and hence requires d — 1 edges. For involu- 
tions with s fixed points, there are (d — s)/2 indepen- 
dent transpositions in its cycle decomposition, and each 
transposition is represented by an edge in the graph. As 
a consequence, q, r and s have to satisfy the following 
condition: 

r{q 2 +q + l- s)/2 = q 2 +q. (32) 

More generally, we define "projective isospectral data" as 
triples (P, {0W},r), where P is a finite projective space 
of dimension at least 2, and {9^'} a set of r nontrivial 
involutions of P, satisfying the following equation 

r(|P|-Fix(0)) = 2(|P|-l), (33) 

for some natural number r > 3. Here Fix(#) = Fix(#W) 
is a constant number of fixed points of P under each 
9^ l \ and |P| is the number of points of P. 

One can now generate all possible pairs of isospectral bil- 
liards whose transplantation matrix is the incidence ma- 
trix of a PG(2, q), with r and q restricted by the previous 
analysis. 

Using the classification of involutions for dimension 2, 
we examine the various cases. 

Let q be even and not a SQUARE. Then any invo- 
lution is an elation and therefore has q + 1 fixed points. 
Therefore, q and r are constrained by the relation 

rq 2 /2 = q 2 +q. (34) 

The only integer solution with r > 3 and q > 2 is 
(r = 3, q = 2). These isospectral billiards correspond 
to the Fano plane PG(2,2) and will be made of d = 7 
copies of a base triangle. 

Let q be odd and not A SQUARE. Then any involu- 
tion is a homology and therefore has q + 2 fixed points. 
Therefore, q and r are constrained by the relation 

r(q 2 -l)/2 = q 2 +q. (35) 

The only integer solution with r > 3 and q > 2 is 
(r = 3,q = 3). These isospectral billiards correspond to 
PG(2,3) and will be made of d — 13 copies of a base 
triangle. 

Let q — p 2 be a SQUARE. Then any involution fixes 
all points in a Baer subplane PG(2,p) and therefore has 
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p 2 +p+ 1 fixed points. Therefore, p and r are constrained 
by the relation 

r{p 4 -p)/2=p i +p 2 . (36) 
There is no integer solution with r > 3 and q > 2. 

Closed circuits. One could also look for isospectral 
billiards with closed circuits: this will require the base 
tile to have a shape such that the closed circuit does not 
make the copies of the tiles come on top of each other 
when unfolded. If we allow just one closed circuit in the 
graph describing the isospectral billiards, then there are 
d edges in the graph instead of d — 1 and the equation for 
p and r becomes 



r(p i -p)/2=p 4 +p 2 + l, 



(37) 



which has the only integer solution (r = 3,p = 2). These 
isospectral billiards correspond to PG(2,4) and will be 
made of d — 21 copies of a base triangle. 

To summarize, we have the following: 

• The Fano plane PG(2,2) provides three pairs 
(made of seven tiles). 

• PG(2,3) provides nine pairs (made of 13 tiles). 

• PG(2,4) provides one pair (made of 21 tiles). 

It turns out that the pairs o btained in s uch a way 
are exactly those obtai ned by iBuser et~al. (1994) and 
Okada and Shudol (|200l[) . 



Now consider the space PG(3, 2), which contains 15 
points. The collineation group of PG(3, 2) is the group 

PGL(4,2) = PrL(4,2) = GL(4,2). (38) 

Generating all possible graphs from the 316 involutions, 
one obtains four pairs of isospectral billiards with 15 
triangular tiles, which c omple tes th e list of all pairs 
found by IBuser et~al. (1994) and lOkada and Shudol 
(|200lh . This list can be found in Appendix \M 



For projective sp aces of d imensi on 2, we thus have the 
following result (|Giraudl . l2005al ). Let P = PG(2,g) 
be the two-dimensional projective space over the finite 
field Fq, and suppose there exists projective isospectral 
data (P,{6>^},r). If q is not a square, then (r,q) € 
{(3, 2), (3, 3)}. If q is a square, then there are no integer 
solutions of Eq. (153")) . 

So the method introduced by iGiraud ( 2005a ) explicitly 
gives the transplantation matrix T for all these pairs — 
each one is the incidence matrix of some finite projec- 
tive space, and the transplantation matrix provides the 
mapping between eigenfunctions of both billiards. The 
inverse mapping is given by 

(l/q n -')(T T -(\/k)X). 



2. Generalized isospectral data 



(39) 



Thasl ( 2006a ) obtained the next generalization for any 
dimension n > 2. It turns out that all possible candidates 
PG(n, q) other than the ones already obtained are ruled 
out by the following results. 



Theorem IV.4 (|Thasl (f2006a^ Let P = PG(n,?) be 

the n- dimensional projective space over the finite field 
¥ q , and suppose there exists projective isospectral data 
(P, {9^}, r). Then q cannot be a square. If q is not a 
square, then (r,n,q) € {(3, 2, 2), (3, n, 3)}, where in the 
the case (r, n, q) — (3, n, 3) each flW fixes pointwise a hy- 
perplane, and also a point not in that hyperplane. How- 
ever, this class of solutions only generates planar isospec- 
tral pairs if n — 2 . 

Call a triple (P, {9^}, r), where P is a finite projective 
space of dimension at least 2, and {9^} a set of r non- 
trivial involutory automorphisms of P, satisfying 



|P|)-£Fix(0«) = 2(|P|-l) 



(40) 



for some natural number r > 3, "generalized projective 
isospectral data" . 

These data were completely classified in (|Thasl . l2006bh . 



Theorem IV.5 (|Thasl (|2006bh ) Let P = PG(l,q) be 

the I -dimensional projective space over the finite field ¥ q , 
I > 2, and suppose there exists generalized projective 
isospectral data (P,{9^},r) which yields isospectral bil- 
liards. Then either I = 2, the 9^ l > fix the same number of 
points o/P, and the solutions are as previously described, 
or I = 3, r = 3 and q = 2, and again the examples are as 
before. 



3. The operator group 

The same kind of results can be formulated at a more 
abstract level. Suppose D is a Euclidean Tl-domain 
on d base triangles, and let fi £ {1,2,3}, be 

the corresponding permutation d x ri-matrices. De- 
fine again involutions on a set X of d letters 
Ai, A 2 , . . . , Arf (corresponding to the base triangles) 
as follows: 9 ( -^(A l ) = Aj if M^f = 1 and i ^ j. In 
the other cases, A, is mapped onto itself. Then clearly, 
(0vW || /i g {1,2,3}) is a transitive permutation group 
on X, which we call the operator group of D. 



Suppose that (Di,D2) is a pair of non-congruent planar 
isospectral domains constructed from unfolding an r-gon, 
r > 3, d < oo times. Since Di are constructed by unfold- 
ing an r-gon, we can associate r involutions 9\ to Di, 



J = 1,2,. 



and i = 1,2. Define the operator groups 



Now suppose that 



G 1 =PSL(n,cy) = G 2 , 



(41) 



(42) 



with q a prime power and n > 2 a natural number. The 
natural geometry on which PSL(n, q) acts (faithfully) is 
the (n— l)-dimensional projective space PG(n— 1, q) over 
the finite field ¥ q . It should be mentioned that PSL(n, q) 
acts transitively on the points of PG(n— 1, q). So we can 
see the involutions 6^ for fixed i G {1,2} as automor- 
phisms of PG(n — 1, q) that generate PSL(n, q). 
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This implies (by nontrivial means) that for fixed i € {1,2} 
the triple 



(PG(n-l,q),{6 ( f ) },r 



(43) 



yields generalized projective isospectral data. The- 
orem IIV.5I implies that (n, q) is contained in 
{ (3, 2), (3, 3), (4, 2), (3, 4)} if n> 3. 

Now suppose that n = 2. We have to solve the equation 

r 

r|PG(l,?)| - VjFix^P) = 2(|PG(l,g)| - 1), (44) 

3=1 

for fixed i € {1,2}, where Fix(# J -' J ' ) ) is the number of fixed 

points in PG(1,?) of 0^ ■ Since |PG(l,g)| = q + 1 and 
since a nontrivial element of PSL(2,g) fixes at most 2 
points of PG(1, q), an easy calculation leads to a contra- 
diction if q > 3. 

Now let q = 2. Then PSL(2, 2) contains precisely 3 
involutions, and they each fix precisely one point of 
PG(1,2). A numerical contradiction follows. ■ 

Thus, the only possible examples of isospectral bil- 
liards that can be construct ed from the third f amily 
of finite simple gro ups (see (Conway et all 1985|)) are 



those obtained in dBuser et all 11994 iGiraudl l2005bl: 



Okada and Shudol . 120011) . They are listed in Appendix 



A. 



V. SEMICLASSICAL INVESTIGATION OF 
ISOSPECTRAL BILLIARDS 

The existence of isospectral pairs proves that the knowl- 
edge of the infinite set of eigenenergies of a billiard does 
not suffice to uniquely determine the shape of its bound- 
ary. A natural question arises: if the set of eigenvalues is 
not sufficient to distinguish between two isospectral bil- 
liards, then which additional quantity would suffice to 
uniquely specify which is which? A parallel issue is to 
identify what kind of geometric information on the sys- 
tem one can extract from the spectrum. This type of 
inverse problem occurs in many fields of physics, from 
lasing cavities to stellar oscillations. 
It is well known that classical mechanics can be seen as 
a limit of quantum mechanics when Planck's constant, 
seen as a parameter, goes to zero. It is therefore natural 
that, for small enough values of this parameter, classical 
characteristics of quantum systems begin to emerge. If 
one considers an electron in a box, one can construct a 
certain linear combination of stationary wave functions 
that describes its probability density distribution at each 
point of the box. At the classical limit, this probabil- 
ity distribution gets mainly concentrated on classically 
authorized trajectories. The quantum system thus some- 
how "knows" about classical trajectories of the underly- 
ing classical system. As shown in this section, the semi- 
classical approach provides a constructive way to retrieve 
geometric information on the system. 
More formally, the time-dependent propagator of the 
Schrodinger equation can be expressed as a Feynman 
path integral, which is a sum over all continuous paths 
going from the initial to the final point. Using a 



stationary phase approximation, IVan Vleckl (|l928l ) ob- 
tained a formula expressing the propagator (or, more pre- 
cisely, its discretized version) in the semiclassical limit 
as a sum over all clas sical trajectories of the system. 
iBalian and Bloch ( 1974 ) showed that the density of states 
can be written as a sum over closed trajectories of the 
classical system. Using a stationary phase approxima- 
tion technique, the semiclassical Green function can be 
similarly expressed as a sum over all classical trajecto- 
ries. This led to the Gutzwiller trace formula for chaotic 



systems (see (jGutzwillerl . 119911) and references therein) 
or the Berry- Tabor trac e formula for integrable systems 
( Berry and Tabor , 19761) . These trace formulae relate the 
quantum spectrum to classical features of the system. 
While the leading terms of the mean spectral density pro- 
vide geometric information about global quantities of the 
system, such as the area or perimeter, the trace formulae 
contain information about classical trajectories. Correc- 
tions to these trace formulae account for the presence of 
other classical trajectories, such as diffractive orbits. 
As mentioned in section III. At the transplantation proof 
of isospcctrality shows that pairs displaying any kind 
of classical dynamics can be constructed, from (pseudo- 
)integrability to chaos. One might ask whether the spec- 
trum of a billiard uniquely determines its length spec- 
trum. As we will see, the transplantation method pro- 
vides an answer to this question. However, in the pseudo- 
integrable case where diffractive contributions to the 
trace formula can be handled, it turns out that transplan- 
tation properties of diffractive orbits are different from 
those of periodic orbits. 

In this section we first introduce some tools relevant to 
semiclassical quantization and then review in more detail 
various classical and quantum properties of isospectral 
pairs that have been studied in the literature, either for 
generic isospectral billiards, or for particular examples 
such as the celebrated example of Fig. [TJ 



A. Mean density of eigenvalues 

The problem of calculating the eigenvalue distribution for 
a given domain B (sometimes called Weyl's problem) is 
dealt with starting from the density of energy levels 



d(E) =J2 S ( E - E n 



(45) 



where S is the Dirac delta function and the sum runs 
over all eigenvalues of the system. The counting function 
is the integrated version of the eigenvalue distribution: 



N{E) = ©(£ 'En), 



(46) 



where O is the Heaviside step function. Statistical func- 
tions of the energy can be studied by proper smoothing 
of the delta functions in (|4"5"j) . The mean of a function 
/ of the energy is defined by its convolution with a test- 
function £: 



f(e)£(E - e)de. 



(47) 



The test-function £ is taken to be centered at 0, normal- 
ized to 1 and have an important weight only around the 
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origin, with a width AE large compared to the mean level 
spacing but small compared to E. 

Isospectral billiards share by definition the same counting 
function J\f(E). Let us study the mean behavior of JV(E). 
Suppose the Hamiltonian of an iV-dimensional system is 
of the form 



certain relation between the angles cdi of the first billiard 
and the angles o! i of the first billiard must hold, namely 



first billiard 



E 

second billiard 



7T 



(54) 



H(q,p)=p 2 /2m + V(q). 



(48) 



The "Thomas-Fermi approximation" consists in making 
the assumption that each quantum state is associated 
with a volume (2irh) N in phase space. The mean value 
of N(E) is given by 



1 / m \ N / 2 



r(f + i)V2 7 rfi2 / Jv[q)<E 



(49) 

[E - V{q)] N ' 2 d q 



after integration over p. In the case where we describe 
the movement in an TV-dimensional domain of volume V 
we get 



77(E) 



V 



) E N ' 2 



r(j\r/2 + i) \2ith 2 ) 



(50) 



which is the first term in a series expansion of 77(E), 
called the Weyl expansion. In particular two isospectral 
iV-dimensional domains necessarily have the same vol- 
ume. 

For two-dimensional billiards and under our conventions 
on units, this first term of Weyl expansion reads 



77(E) * ±E, 



(51) 



where A is the area of the billiard. This means that 
a necessary condition for isospectrality is that the bil- 
liards have the same area. The asymptotic expan- 
sion of the Laplace transform o f the density of states 
(jStewardson and Waechten . Il97ll ) allows us to derive the 
follow ing terms in the Weyl expansion (jBaltes and Hill . 
197G). The expansion is given by 



AC 

77(E) ^ -p-E =f j- 

4-7T 47T 



E + K, 



(52) 



where A and C are the area and the perimeter of the 
billiard, respectively. The sign before C is (-) for Dirich- 
let boundary conditions and (+) for Neumann boundary 
conditions. The constant K, depends on the geometry of 
the boundary. For boundaries with smooth arcs of length 
7, and corners of angle < ay < 27r it reads 



24 V a t / ^— ' / Ott ' 



«(0. 

2tt 



(53) 



where n(l) is the curvature measured along the arc. 
The Weyl expansion ([52"T) shows that if two billiards have 
the same spectrum, then they necessarily have the same 
area and the same perimeter. Furthermore, a certain 
combination of the properties of their angles and curva- 
tures must be the same. In the case of polygonal isospec- 
tral billiards, such as those given in the examples in Ap- 
pendix El the fact that JC must be the same entails that a 



B. Periodic orbits 

The previous section gives necessary relations that must 
hold between two isospectral billiards, in particular the 
fact that they must have the same area and perimeter. 
These relations were based on the fact that the mean 
density of quantum eigenvalues (or the mean counting 
function) could be related to classical features of the bil- 
liards. In fact deeper relations exist between the quantum 
properties of a billiard and its classical features. These 
relations are expressed through "trace formulas", which 
express the density of energy levels as a sum over classical 
trajectories, in the semiclassical approximation. Semi- 
classical methods are based on the fact that the classi- 
cal limit of quantum mechanics is obtained for h — > 
in the path integral expressing the propagator. The ex- 
pansion of this path integral in powers of h allows us to 
calculate the sequence of quantum corrections to classi- 
cal theory. The semiclassical approximation keeps in this 
expansion only the lowest-order term in h. Corrections 
to this approximation correspond to taking into account 
higher-order terms. In this section we recall the main 
steps leading to a trace formula for billiards, and apply 
it to isospectrality. 



1. Green function 

The propagator of the system is defined as the conditional 
probability amplitude K(qf,tf,qi,ti) for the particle to 
be at point qj at time tt, if it was at point qi at time ti. 
The propagator is the only solution of the Schrodinger 
equation that satisfies the condition 



lim K(qf,tf,qi,ti) = 6(qf-qi). 



(55) 



One can then show that the propagator can be written 
as a Feynman integral 



K(Qf,tf,qi,U 



Vq(t)e 



jr f dtL{q,q,t) 



(56) 



where the sum runs over all possible trajectories going 
from (qi,U) to (qf,tf) and L is the Lagrangian. The no- 
tation ([551) has to be understood as the limit as n goes 
to infinity of a discrete sum over all n step paths go- 
ing from (qi,ti) to (qf,tf): the integral (|56|) runs over 
all continuous, but not necessarily derivable, paths. One 
immediately sees that the classical limit of quantum me- 
chanics corresponds to letting the constant h go to 0: 
the main contributions to the probability K then cor- 
respond to stationary poin ts of the action J dtL(q,q,t) 
( Feynman and Hibbsl . 1965 ). 

The advanced Green function is the Fourier transform of 
the propagator, which is defined by 



G{q f , qi ;E) 



ih 



dt K(qf,t; qi ,0) e lEt ' h . (57) 
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It is a solution of 

(-H + E)G(q f ,q i ;E)=S(q f -q i ). (58) 

The action along a trajectory can be denned as the inte- 
gral of the momentum 



rif 

S(q/,qi',E) = pdq, 

Jqi 



(59) 



and the Green function as 



1 



G(q f , qi ;E) = ^J Vq(t)e^^ E \ (60) 

where the path integral now runs over all continuous 
paths going from qt to qf at a given energy E. 
In many cases Eq. (|58|) allows us to calculate the Green 
function. In the case of free motion in Euclidean space, 
the Hamiltonian reduces to the Laplacian (up to a sign), 
and the Green function is solution of 



(A, + E)G(q f , qi ; E) = 6{q f - 



(61) 



where the qf index recalls that the derivatives of the 
Laplacian are applied on variable qf. In two dimensions, 
the Green function reads 



1 



G{q f ,q i ;E) = j,H£\k\q f -q i 



with k = V E and Hq the Hankel function of the first 
kind. 



(62) 



2. Semiclassical Green function 

The expression for the Green function G(qf,qi;E) 
is a sum over all continuous paths joining qi to qf at 
energy E. The semiclassical approximation consists in 
keeping only the lowest-order term in the h expansion. 
This term is given by the stationary phase approximation. 
The only paths contributing to integral (150)) are paths 
for which the action S reaches a stationary value, that 
is, paths that correspond to classical trajectories. The 
semiclassical Green function can thus be expressed as a 
sum, over all classical trajectories. Each term in the sum 
is an exponential whose phase is given by the classical 
action integrated along the trajectory. The prefactor is 
obtained by the stationary-phase approximation around 
the classical trajectory. 

Choosing a coordinate system (q\\,q±) such that gy is the 
coordinate along the trajectory and q± is the coordinates 
perpendicular to the trajectory, one obtains the semiclas- 
sical Green f unctio n as a sum over all classical trajectories 
(|Gutzwillerl . ll99l[ ): 



G sc ( qf , qt ;E) = J2 



2n 



(2inhY N+1 '>/ 2 



(63) 



i I d 2 s 

— dct -- 



1 1/2 



* II 9/11 



dqf ± dq t 



1 7T 

x exp \ -S(qf,qr,E) - 



where N is the space dimension. The phase /i is called 
the Maslov index of the trajectory. In two dimensions 



for hard wall reflections, each reflection of the classical 
orbit yields a contribution [i = 2 for Dirichlet boundary 
conditions and /i = for Neumann or periodic boundary 
conditions. 



3. Semiclassical density of eigenvalues 

We defined the Green function of a quantum system by 
Eq. ([57]) . It will be more useful to express the Green 
function as a sum over eigenvalues and eigenfunctions of 
the Hamiltonian. It can be verified that formally 



G(q f ,q i ;E) = Y J 



^»(gi)^n(g/) 

E — E n 



(64) 



where ^ denotes the complex conjugate of is a solution 
of Eq. ([58l) . In order to give a mathematically correct 
meaning to this expression, we use the advanced Green 
function 



G+(qf,qi;E) = G{q f ,qi;E + ie). 



(65) 



The words "Green function" will always implicitly refer 
to the limit of the advanced Green function for e — > 0. 
The density of energy levels (l4"5j) can be related to the 
Green function by 



d(E) = -~flm G{q,q;E) 



dq. 



To prove this, we use the fact that for e — > 0, 

lim = P it:8(x) 

e->o x + ie x 



(66) 



(67) 



(P denotes the principal value and 8 is the Dirac delta 
function), and that, since H is Hermitian, its eigenvectors 
verify J vE'm^E'n = S mn . The Green function G(q' , g; E) di- 
verges for q' — >■ q but not its imaginary part. The expres- 
sion Im G(q, q; E) has to be understood as the imaginary 
part of G(q', q) taken at the limit q' — > q. Thanks to this 
relation, the density of states can be expressed as the 
trace of the Green function. Equation (|66|) is the start- 
ing point of trace formulae. Note that, if the density of 
states (|45[) is regularized as a sum of Lorentzians 



d e (E) = ~Y / 

7T ^ ^ 



1 



(e - E n y + e 2 



one gets 



d e (E) 



Im G(q, q;E + ie) dq. 



(69) 



Equation (|66| must therefore be understood as the limit, 
as e — > 0, of each member of Eq. (|69|) . However, the 
density of states is usually calculated from the Green 
function by first evaluating the integral for q = q' (the 
"trace" of the Green function), and then taking the 
imaginary part. This can be made rigorous, by multi- 
plying the Green function by some factor making the 
integ ral convergent in the limit q = q' (jBalian and Blochl . 
1974 ). 
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The semiclassical density of states is then obtained by 
use of Eq. (frJ3"j) with g^=g , /. The density of states in the 
semiclassical approximation is then the sum of a "smooth 
part" and an oscillating term that is a superposition of 
plane waves, 



d s - c -{E) = d(E) + d osc (E). 



(70) 



The term d is obtained from the first term (|5Tj) of the 
Weyl expansion. It gives a mean density of states equal 
to 



7 A 

a = — . 

47T 



(71) 



The oscillating term reads 



d osc (E) 



( 2 ^) 3/2 ptrJdet(M» 



'(t-^t ). 



-c.c. 



(72) 

The Gutzwiller trace formula (|72p is a sum over all prim- 
itive periodic orbits (pp), repeated n times. Each prim- 
itive periodic orbit has a certain action S p , period T p , 
monodromy matrix M p , and Maslov index v v (taking into 
account additional phases owing to integration). The 
identity matrix is denoted by I, and c.c. denotes the 
complex conjugated. 

In the case of integrable and pseudo-integrable systems 
(such as the isospectral pair of Fig. [T]), periodic or- 
bits are no longer isolated but appear within families of 
parallel trajectories having the same length ("cylinders 
of periodic orbits"). The Gutzwiller trace formula no 
longer applies. Pseudo-integrable billiards are both non- 
integrable and non-chaotic, and their classical character- 
istics are intermediate between those of integrable and 
those of chaotic billiards. Classical trajectories appear 
within families of parallel trajectories of same length, 
but nevertheless the equations of motion are not exactly 
solvable because of the presence of diffraction corners. 
Berry and Tabor derived a trace formula for mul- 

tidimensional integrable systems that can be adapted to 
polygonal billiards. In the case of a two-dimensional 
polygonal billiard, the trace formula becomes 



d osc (E) 



pp 



A 

71=1 



9^ 



-3i7r/4— irtVpp-ir /2 



+ c.c, (73) 



where A p is the area occupied by the cylinder of periodic 
orbits labeled by p. Equation ((73)) gives us a strong 
relationship between periodic orbits of billiards having 
the same spectrum. The trace formulas must be the 
same, and one might think that the equality of the sums 
over periodic orbits can be achieved only if the periodic 
orbits are identical in the two billiards. 

It turns out that this is true. It can be proved fairly eas- 
ily that two transplantable isospectral domains have the 
same length spectrum (i. e. both domains have pe riodic 
orbits of the same length) (jOkada and Shudol . l200ll ). The 
proof is given in section lV.I.ll Here we illustrate this fact 
on a simple example. Consider the billiards of Fig. [71 It 
is possible to encode any trajectory drawn on the billiard 
(provided it does not pass through vertices) by symbolic 
dynamics. Consider a trajectory 77j drawn on the first 



billiard, going from tile i to tile j. Recall that the way 
the building blocks are glued together (or, equivalently, 
the coloring of the associated graph) can be described 
by matrices M^\N^\ 1 < /i < 3, as introduced in 
section [II] With the trajectory %j one can associate a 
"word" (ai, a%, . . . , a n ) describing the sequence of edges 
crossed by the trajectory on its way. To this trajectory 
we then associate the matrix M = ■ Note that 

there exists a trajectory between tiles i and j if and only 
if Mij = 1. We then define N = n^ (ai) - The trans- 
plantation between the two billiards can be described by 
some matrix T such that TM^ = N^T for 1 < /j < 3. 
These commutation relations imply that TM = NT also 
holds. In particular, if k is a tile of the second billiard 
such that Tki — 1 and k' a tile of the second billiard such 
that T k ij = 1, then {TM) kj = 1 = (NT) kj , which implies 
that Nkk' = 1. This is exactly equivalent to saying that 
the trajectory 7y can be drawn on the second billiard 
between tiles k and k'. Fig. [TU] shows two pencils of pe- 
riodic orbits on each billiard. One can check that these 
two pencils appear with the same length and the same 
width in both billiards. 





FIG. 10 Periodic and diffractive orbits in the unfolded pair of 
Fig-El 



C. Diffractive orbits 

The semiclassical trace formula (|72p . which is expressed 
in terms of classical periodic orbits, is only a leading- 
order approximation for small values of h. Higher- 
order corrections to this formula take into account con- 
tributions from diffractive orbits: creep ing trajecto- 
ries, trajectories be t ween s cattering points (IKellerl . Il962 ; 



Pavloff and Schmitl . Il995t IVattav et all . |l29jj), and or 



bits almost tangent to a concave section of the boundary 
(jPrimack et aZ.I . I1997T) . In the case of polygonal billiards, 
the semiclassical trace formula (l73l) has to be corrected to 



take into account scattering trajectories, that is, classical 
trajectories going from one scattering point to another, 
or combinations thereof. 
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As in the case of periodic orbits, one might believe that 
the equality of densities d(E) for isospectral billiards 
must translate to an equality of diffractive orbits. Sur- 
prisingly, this is not the case, as we now show. Again, we 
concentrate on the simple example of polygonal isospec- 
tral billiards. 

In th e case of polygonal billiards, lHannav and Thairi 
( 2003T ) were able to derive an exact expansion for the 
Green function, as a sum over all scattering trajectories. 
The exact Green function between a point a and a point 
b reads 

G(a,b) = Y,j^ £ 1^ j_ oo dsids2 - dSn 



tfe 



x [kR(si,s 2 ,...,s n )} 



n 

k=l 



2tt 



{ lk M k + k + is k ) 2 - 7T 2 ' 



(74) 



where 



R 2 ( Sl ,S2,...,s n ) = (75) 

(r + m si + r 2 e Sl+S2 H h r n e Sl+S2+ - +s ") 

x (r + r x er Sl + r 2 e^ Sl ^ S2 H h r n e~ Sl ~ S2 s ") . 

The Green function appears as a sum over paths made 
of n + 1 straight lines of length r^, < i < n. The first 
line goes from point a to a diffracting corner, then there 
are n scattering trajectories going from one diffracting 
corner to another, and finally a trajectory going from 
one diffracting corner to point b. The diffraction angles 
are M k ^ k + 9 k , 1 < k < n, with j k the measure of the 
angle at the diffracting corner and M k the number of 
times the path winds around the diffracting corner (thus, 

o<e k < ^). 

Giraud (2004) showed, using the expansion (|74[) of the 
Green function, that isospectral domains can be distin- 
guished by the fact that in general the lengths of their 
diffractive orbits differ. This can be illustrated in the 
case of the billiard with rectangular base tile unfolded to 
a translation surface (Fig. [7]). If the sides of the base 
tiles are incommensurate, then there cannot be diffrac- 
tive orbits of the same length as a given diffractive orbit 
but in the same direction in the plane. For instance, for 
the dashed diffractive orbit drawn in the second billiard 
of Fig. [TU1 orbits starting from a diffractive corner of the 
first billiard in the same direction never reach another 
diffractive corner. This means that the dashed orbit has 
no partner in the first billiard. 

The connection between the energy spectrum and the 
length spectrum through the trace formula indicates, 
however, that these discrepancies between diffractive or- 
bits must be compensated in a certain way. This com- 
pensation can be understood by analyzing the formula of 
Hannay and Thain ([74)1 . In fact each contribution to the 
Green function in Eq. (|74|) has to be understood as an in- 
finite sum over all windings around vertices (see Fig.llip. 
Here by vertices we mean the four corners and the two 
points at the middle of the horizontal sides of each of the 
seven rectangular tiles in Fig. [7] If there is a diffracting 
corner (as is the case for instance at the bottom right 
corner of tile 7 in the second billiard of Fig. [7|) then there 
is a non-zero contribution, while if there is no scatterer 



=lim 

e-0 




FIG. 11 A contribution to the Green function in the case of 
forward diffraction. If the orbit goes through a vertex, the 
term in Eq. (|74|l should be interpreted as the limit for e — s- 
of an infinite number of trajectories. If there is no vertex only 
the straight path contribution remains, the other (winding) 
terms add up to zero. 



(e.g. at the bottom left corner of tile 7 in the second 
billiard) the series of diffractive terms adds up to zero, 



00 

E 



2tt 



(2nM k 



is k ) 



= 0. 



(76) 



As a consequence, a diffractive contribution to the Green 
function, going from a point a to a point b through pos- 
sibly several vertices, has to be understood as a sum 
of trajectories winding around both scattering and non- 
scattering vertices (see Fig. [IS]) . Now each of these new 
"fictitious" trajectories avoids vertices (since they wind 
around). The reasoning we used in section IY.BI applies: 
With any trajectory one can associate a matrix M de- 
scribing the edges crossed by the trajectory. The matrix 
N corresponding to M in the other billiard is such that 
TM = NT, and a partner of the diffractive orbit can be 
found between tiles i and j such that Nij = 1. 
Thus, even though diffractive orbit lengths might differ, 
each "expanded" diffractive orbit indeed has a partner of 
the same length. 




terms of 
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FIG. 12 A contribution to the Green function in the case of 
forward diffraction. The filled circles are scattering vertices; 
the empty ones are non-scattering vertices. 



D. Green function 

This relation between diffractive orbits translates to a 
relation between Green functions of the two domains 



(Giraud, 2004). The matrices M,N introduced in sec- 
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tion IV.BI verify the property 



(77) 



which can be proved using the commutation relation (|7]l 
and the fact that (T 2 )^ = 1 + 2<5y. Thus, in the expan- 
sion (174)) of the Green function between a point in tile i 
and a point in tile j in the first (second) billiard, each tra- 
jectory appears with a weight My {Nj)- But according 
to Eq. J77| we have 



u 2 



o o ' 



(78) 



Therefore from Eq. (f74")) and identity ((75)) one can infer a 
relation between Green functions, namely, 

G^(a,t;b,j) = ^T u >T jj .G^\a ) i';bJ , )~G^(a-,b), 

(79) 

where G^\a;b) is the Green function on the base tile. 
This relation between Green functions like the relations 
between periodic orbits or diffractive orbits are all con- 
sequences of the transplantation property which is the 
fundamental feature of all known examples of isospectral 
billiards. 



E. Scattering poles of the exterior Neumann problem 



In section fV.CI we considered the particular case of polyg- 
onal isospectral billiards, for which it is possible to ex- 
press the exact Green function as an infinite expansion 
given by Eq. ([71]) . In a more general setting, it is also 
possible to express the Green function of the billiard 
with Dirichlet boundary conditions as an infinite sum 
taking into accou n t all p ossible reflections on obstacles. 
Balian and Blochl (|1974l ) gave a general method, called 
"multiple reflection expansion", which gives the Green 
function in terms of the free Green function Go . Applied 
to a two-dimensional billiard, this expansion is given by 



G(q,q';E)=Go(q,q';E) 



(80) 



-2/ ds Go(q,s;E)d s G (s, q';E) 

JOB 

+ (-2) 2 / ds ds' d s G (q,s;E) 
JdB 



x d s ,G (s,s';E)G (s',q';E) + --- , 

where s and s' are points along the boundary, and d x 
denotes the derivative along an outward vector normal 
to the boundary at point x. The first term Go(q,q'; E) 
on the right-hand side of Eq. (|H0"| corresponds to direct 
(free) propagation from q to q' ', the first integral corre- 
sponds to trajectories from q to q' with one reflection on 
the boundary at point s, and so on. We introduce the 
kernel KE{q,q') = —2d q 'Go(q,q';E), which is a contin- 
uous infinite-dimensional operator defined on dB x dB 
(dB is the boundary of the billiard). One can express 
Eq. ([H as 



G(q,q';E) = G (q,q';E) (81) 

OO „ 

-2 V / ds ds' G (q, s; E)K E (s, s')d s ,G (s', q'; E). 

„-nJdB 



Formal performance of the sum over n yields the infinite- 
dimensional operator (I — K pi) -1 , where I i s the identity 



operator. Fredholm theory ( Smithieslll962t ) showed that 
for sufficiently "nice" billiards the operator (I — K E ) 
well-defined and can be expressed as 



is 



(i - K E y 



D{EY 



(82) 



where D(E) is the Fredholm determinant det(I — K E ), 
and Ne is an infinite-dimensional operator defined on 
dB x dB. The Fredholm determinant admits an expan- 
sion 



D(E) = J2 D n (E), 



(83) 



n=0 



with Dq(E) — 1, and for n > 1 

(-iy 



D n (E) = 



dqi 



OB 



dB 



dq n K(q,q) (84) 



(q is the vector (<?i, q 2l . . . , q n ))- We have introduced the 
determinant 



K(q,q') = 



K E {qi,q' x ) K E (q 1 ,q' 2 ) ■■■ K E {q\,(f n ) 
K E (q2,q'i) K E (q 2 ,q' 2 ) ■■■ K E (q 2 ,q' n ) 

K E {q n ,q'i) K E (q n ,q' 2 ) ■■■ K E {q n ,q' n 



(85) 

The operator Ne is defined on dB x dB by its expansion 



N n = J2Dk(E)K£- k . 

k=0 



(86) 



The Fredholm determinant D(E) appearing in the ex- 
pression of the Green function has the property that it has 
zeros at eigenvalues of the system (jGeorgeot and Prangel 
1995T ). A natural question is whether isospectral billiards 
share the same Fredho lm determina nt. 
It has been shown by iTasaki et all 



(|l997f ) that, for bil- 
liards with C 2 boundary, D(E) can be decomposed into 
an interior and an exterior contribution, namely D(E) = 
D(0)di nt (E)d ex t(E) . The exterior contribution d ex t(E) is 
related to the scattering of a wave on an obstacle having 
the shape of the billiard with Neumann boundary condi- 
tions, i.e. the zeros of its analytic continuation are res- 
onances of the exterior scattering problem. The interior 
contribution reads 



L 2 E 



n=l 



n(i-# 



E n 



,E/E n 



(87) 

where A and C are the area and the perimeter of the 
billiard, respectively, and 7 is a constant depending on 
the geometry of the billiard. The zeros of d- ul t(E) are 
thus the eigenenergies of the interior Dirichlet problem. 
Obviously, isospectral billiards share the same interior 
part di nt (E). But the exterior part depends on the 
shape of the billiard. In particular solutions of the 
exterior Neumann scattering problem may differ between 
two isospectral bil liards. Therefore a conclusion of 
(jTasaki et all , Il997l) is that isospectral pairs might be 
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distinguished by measuring the sound scattered by them. 

To check t his property, numeric al investigations were per- 
formed bv lOkada et al. ( 2005a ). In fact, Fredholm theory 
applies only for billiards with a smooth boundary, which 
is not the case for any of the known examples of isospec- 
tral pairs. For billiards with a piecewise smooth bound- 
ary, it is however possible to approximate the Fredholm 
determinant D(E) by a discretized version D m (E), de- 
pending on the number m of points taken on the bound- 
ary of the billiard, which converges to D(E) for large 
m. This convergence fails for boundaries with corners: 
All D m ( E) tend to 0. Nevertheless, for domains with 



corners (Okada et al 



2005bl ) showed that it is possible 
to define a regularized version of D m (E) that converges 
to D(E) /D(0) = d in t ( E)d ext (E). Using this regularized 
version, Okad a et all (|2005al ) computed numerically ze- 
ros of the regularized Fredholm determinant for various 
pairs of isospectral billiards. It was observed that zeros 
of the determinant close to the real axis coincide, as they 
should since they are eigenvalues of the interior prob- 
lem. On the other hand, complex zeros (remote from the 
real axis), which correspond to resonances of the exterior 
Neumann problem, are shown to differ. To quantify this 
discrepancy between the resonances of the two billiards, 
the resonance counting number 

N s (r) = [z £ C; \z\ < r, ~ < arg(z) < s} (88) 



was studied by Okad a et all (|2005al ). The best fit 
Ns(r) = Cs,rT 2 , computed over the range r £ [0,-R], 
yields noticeably different values of Cs,r for each billiard. 
This clearly shows that isospectral pairs can indeed be 
distinguished by resonances of scattering waves. 



F. Eigenfunctions 

1. Triangular states 

In general, analytical solutions to the Helmholtz equation 
(A + E)^ = with Dirichlet boundary conditions cannot 
be found. However, it is possible to construct particular 
solutions of this equation provided solutions are known on 
elementary subdomains. This is, for instance, the case if 
the subdomains have the shape of a half-square (billiards 
of Fig. [TJ, or a rectangle (billiards of Fig. 
We take the example of the two billiards in Fig. [TJ Each 
billiard is made of seven triangular (half-square) tiles. 
Eigenfunctions for a dx d-square with Dirichlet boundary 
conditions are of the form 



d 2 



sin ( 



mnx\ . fniry\ . „. 



with eigenvalues E m ^ n = 7r 2 (m 2 + n 2 )/d 2 , m,n > 1. 
Eigenfunctions for the elementary triangles with Dirichlet 
boundary conditions are obtained from (I89[) by antisym- 
metrization with respect to the diagonal, 
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TABLE I First triangular modes t m ,n- Left: Dirichlet bound- 
ary conditions. Right: Neumann boundary conditions. 



and the corresponding eigenenergies are given by ir 2 {m 2 + 
n 2 )/d 2 , m > n. For the sake of definiteness, we consider 
the two is ospectral pairs o n a Cartesian reference frame, 
following (|Wu et ali\l99^ , as in Fig. [H The functions 
t m ,n turn out to be elementary solutions of the Helmholtz 
equation for both isospectral billiards of Fig. [T5J Indeed 
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FIG. 13 Isospectral billiards divided into smaller regions. 

tm.n vanishes on all lines x = kd, y = kd, y — x + 
2kd and y = —x + 2kd, k £ Z, which are precisely the 
lines on which the boundaries of both billiards lie (in 
the convention of Fig. I13|) . The particular solutions t mn 
are called "triangular states" . An example of the lowest- 
energy triangular state is given in Appendix [X] 
The labels of the lowest-energy triangular states among 
the eigenvalu es E-\ < Ei < . . . of the billiar ds have been 
calculated bv lGottlieb and McManusI (|l998| ). The results 
are displayed in TableQ] Each integer pair (m, n),m> n, 
defines a triangular state i m , n . Obviously, the fact that 
an integer can be represented in more than one way as 
a sum of two squares leads to degeneracies for triangular 
states, and hence for the isospectral pairs of Fig. [TJ 
Note that for Neumann boundary conditions, it can be 
easily checked that the functions 

4 r fi 

1 V d J V d 
'rmty\ fnirxy 
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hr) cos (- 
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(91) 



d J V d 

for < m < n, (m,n) ^= (0,0), have a normal deriva- 
tive that vanishes on all lines x — kd, y — kd, y — 
x + 2kd, and y = —x + 2kd, k £ Z. Therefore u m ^ n 
are solutions of Helmholtz equations for the billiards of 
Fig. [T3] with Neumann boundary conditions. Their label 
among the eigenstates of the b illiards is given in Table [T] 
([Gottlieb and McManusI . Il998t ). 



(90) 



2. Mode-matching method 

The knowledge of these particular triangular states is the 
starting point for the so-called "mode-matching method" . 
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It consists in dividing the billiards into subdomains for 
which solutions of the Helmholtz equation are known ana- 
lytically. Consider for example the left billiard of Fig. [TBI 
It is made of five elementary domains, three triangles 
A, B, E and two squares C, D. For each subdomain, an- 
alytical solutions for the Dirichlet problem are given by 
(translations of) functions (l89l) or (|90l) . We dehne the 
function ip n (x,y) — sin(a n x) sui(b n y) / sin(6 n d), where d 
is the length of the side of the elementary square, and 
we have set a n = nir/d and b n = y/E — a„. If given 
boundary conditions are imposed on the boundaries of 
these subdomains, as in Fig. [TU solutions can be written 
explicitly for these elementary subdomains as superposi- 
tions of functions obtained from translations or reflections 
of tp n - In particular, one can construct functions taking 
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FIG. 14 Elementary regions building the isospectral pairs. 

the value on the plain boundary and sin(a„x) (sin(a„y)) 
on the dashed horizontal (vertical) boundary for each of 
the domains shown in Fig. 1141 Such functions are given 
by 

<$> (A) {x,y) =r/> n (x,d-y)-i/} n {y,d-x) (92) 
®n fay) = i>n{%,y) - ip n (d- y,d- x) 
^n^fav) = i>n(y,x) - ipn(d - x,d - y) 

$l cl) (z,y) = Mv,x) 

^ 2 \x,y)=Mv,d-x) 
^ n D1 \x,y)=^ n {y,d-x) 
$i m Xx,y)=^ n fad-y) 
®nfay) = i>n{x,y) - ip n (y,x). 

The mode-matching method consists in looking for a so- 
lution \P of the Helmholtz equation as a superposition of 
such functions, with amplitudes chosen such that >F and 
its partial derivatives be continuous at each boundary be- 
tween subdomains. At the boundary between elementary 
subdomains, the eigenfunction can be expanded on the 
functions ip n (x) = sin(a n x) as 



where the sum goes from 1 to some trunca- 
tion number N. The eigenfunction $ is en- 
tirely determined by knowledge of the vector 
V = (Ai, . . . ,An, Bi, . . . ,Bff, Ci, . . . , Cjy, Di, . . . , Djy). 
Therefore \& can be written as 

^ A {x,y) = Y,A n ^ n A) {x,y-2d) 

9 B fa y) = A»*i m) fa y-d) + J2 B,M n B2) fa y - d) 

n n 

* c fa y) = Y, B ^n C1) (x-d,y-d) 

Th 

+ Y,Cn^ 2) {x-d,y-d) 

Th 

*ntx,y)=Y J Cn$ ( n D1) (x — 2d,y — d) 
+ J2D n <f { n D2) (x-2d,y-d) 

71 

^ E (x,y) = Y,D 1 M n E \x-2d,y), (94) 

n 

where ^>x is the restriction of the function <3> to the el- 
ementary domain X = A, B, C, D, or E. The function 
^ is indeed an eigenfunction of the billiard if its normal 
derivatives at the boundaries between domains are con- 
tinuous. This latter condition can be written as a sys- 
tem of linear equations that can be cast under the form 
AfV = 0, where M is a AN x 47V-matrix given by 
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(95) 

with U mn = (b n cotb n d)S mn , V mn = (b n / smb n d)S mn , 
W mn = a m a n /(E - c? m - a%), and P mn = {-l) n S mn . The 
matrix M depends on E through b n . Eigenvalues of the 
billiard correspond either to values of E where det M = 
or to V = 0. In the case V — 0, the wave function 
vanishes on the boundaries between the domains, and 
the eigenfunction is a triangular state. If det M = 0, 
Eqs. (|94|) give the corresponding eigenfunction. 



Interestingly, the mode-matching method provides an al- 



ternative proof to isospectrality (|Wu et all 11995). The 
matrix M' corresponding to M for the right billiard of 
Fig. [13] is the AN x 47V-matrix given by 



M' = 



^AB(x,y) 


= ^A n ip n (x) 

n 




(93) 
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It can be easily checked that M and M' are related by 



T = 



M =* TM'T, 
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If \I/ is a solution of the Helmholtz equation for the first 
billiard, it can be written under the form (IM|) with con- 
stants specified by some vector V verifying MV = 0. 
Let be the function defined on the second billiard by 
some constants given by the vector V = TV. Because of 
Eq. (j97l) the vector V verifies M'V = 0, and therefore 
'$>' is a solution of the Helmholtz equation for the second 
billiard. Since the relation between "J and W is linear, 
the cigcncncrgy is the same for both functions, and thus 
the billiards are isospectral. 



G. Eigenvalue statistics 

As explained in section [Til the shape of the elemen- 
tary building block of a pair of isospectral billiards can 
be varied at will provided some conditions are satisfied. 
Thus, examples of chaotic pairs, or pseudo-integrable 
pairs, or even pairs with a fractal boundary can be pro- 
duced. However, the most popular examples of isospec- 
tral billiards, e.g., those of Fig. [TJ are constructed with 
a triangular-shaped base tile. The resulting billiards 
are thus polygonal billiards. Billiards with a polygonal 
boundary can display a whole range of classical behaviors 
from integrability to chaos. Isospectral billiards made of 
tiles whose angles a re rational multiples of tt are pseudo- 
integrable billiards (jRichens and Berrvl 119811 ) . The prop- 
erties of these billiards were mentioned in section Hi. Bl 
In the field of quantum chaos, many works have been 
concerned with a characterization of the statistical prop- 
erties of spectra of billiards. The question of the spectral 
properties displayed by polygonal isospectral billiards at- 
tracted some interest in the literature. Eigenvalue statis- 
tics for th e pair of Fig. [T]have been studied numerically by 
Wu et all (|1995l ). based on the first 598 energy levels. The 
short-range correlations of the spectrum were shown to 
lie between the random matrix statistics of the Gaussian 
orthogonal en semble (GOE) and Poisson statistics (see 
' Ported. Il965l) for a review on the seminal papers, and 



Guhr et al for a recent review on random matrix 



theory). On removal of the 78 triangular states, it was 
observed that the nearest-neighbor level spacing distribu- 
tion function P(s), which characterizes the distribution 
of the spacings between nearest-neighbor energy levels, 
agrees with the nearest-neighbor distribution f or GOE 
matr ices. The spectral rigidity A 3 (L) (see, e.g., (IMehta 



1990) for a rigorous definition) measures the deviation 
of the integrated density of states N(E) (the number of 
eigenvalues smaller than E) from a straight line, on an in- 
terval [E— L/2, E+L/2}. Computation of A 3 showed that 



i t is a lso of GOE type for these billiards. lAurich et al 
(.1997*) calculated the functions E{k,L), which give the 
probability to find k energy levels in an random interval 
of length L ( Aurich and Steined . ll990f) . for isospectral bil- 
liards shaped as in Fig. [TJ again showing a behavior that 
is intermediate between the chaotic and integrable cases. 



H. Nodal domains 

Nodal lines for two-dimensional billiards are one- 
dimensional curves on which eigenfunctions vanish. 
Nodal domains are connected regions of the billiard where 



an eigenfunction has a con stant sign. A theorem by 
Courant and Hilberd (|1953l ) states that the nth eigen- 
function ^„ has at most n nodal domains. The num- 
ber v n of nod al domains in ^„ can be further estimated 



(|Pleiiell . ll956h . We define a rescaled nodal-domain num- 
ber £ n = v n jn G [0, 1]. If ji is the first zero of the Bessel 
function Jo, then lim sup,^^ < (2/ji) 2 . The limit 
distribution of £„ is defined by 



P(0 = lim 
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E-foc Nr 

1 C 



[E] 



E 



(99) 



where I g (E) is the interval [E, E + gE] for some fixed 
g > 0, and Nj is the numbe r of eigenvalues in t he interval 
/. It has been shown by iBlum et al\ (2002) that this 
distribution has universal features. 

For some instances of i sospectral pairs, such a s flat 
tori i n R" with n > 4 ( Gnutzmann et all , 120051 ) (see 
also ( Levitin et al , 20061 )). it was conjectured that two 
isospectral domains produce a different number of nodal 
domains (domains separated by nodal lines where ^ = 0). 
Heuristic argumen ts as well as nume r ical i nvestigations 
were collected by iGnutzmann et al. ( 2005 ) to support 
this conject ure. A recent so l ution of this conjecture can 
be found in (jBruening et all 120081 ). 



I. Isospectrality versus isolength spectrality 

We now consider a related important question. Since 
transplantation is a mapping between the two billiards, 
the classical properties should map onto one another as 
well. Here we investigate the mapping between periodic 
orbits. 



1. Okada and Shudo's result on isolength spectrality 

Let D be a planar domain obtained by unfolding N times 
the same triangular building block B with sides 1,2,3. 
Then the length spectrum is the set of lengths of closed 
trajectories (periodic orbits) of D. Any periodic orbit on 
D can be regarded as a "lift" of a closed trajectory on 
B, because its projection is always a periodic orbit on B. 
(The converse is, of course, not necessarily true.) One 
observes that the number of closed lifts of a given closed 
trajectory on B is counted as 



'(7) = TV {M (t™ ) Af (T" 



.M< 7l) ), 



(100) 



where 7 = Yili (T< ^ {1,2, 3}) denotes the sequence 
representing the order in which a given closed trajec- 
tory on B hits the boundary segments. (The M( 7j )'s 
are adjacency matrices.) Note that such a sequence is 
not uniquely determined by a given closed orbit — the 
number of closed lifts, however, is. So the length spec- 
trum of D is determined by the length spectrum of B and 
by n D (7). Hence, if one considers two domains D and D' 
that are constructed by unfolding the same building block 
as above, it is sufficient to prove that n D {^) = n D (7) 
for all possible sequences 7 in order to deduce "isolength 
spectrality" . 

The following is now obvious. 
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Theorem V.l (jOkada and Shudol (|200lh ) Let D 

and D' be two unfolded domains obtained by N times 
successive reflections of the same building block. If D 
and D' are transplantable, then n D (j) = n D (7) for any 
sequence 7, so D and D' are isolength spectral. ■ 



fP and p respectively, the representations are similar. So 
there exists an invertible matrix T for which 



for any p. Thus D and D' are transplantable. 



(108) 



Let S be a finite set, say, S = {ai, . . . , at} with k G No- 
The free group F = F(S) generated by S is defined as 
follows. Introduce a set S^ 1 := {a J" ,...,Oj } which 
consists of the "inverse symbols" of S. A word with al- 
phabet S (or SUS -1 ) is just a finite sequence of elements 
of S U S . A reduced word is a word in which any se- 
quence consisting of an element of S and its inverse is 
deleted. By definition, F consists of all reduced words 
with alphabet S, together with the empty word. Group 
operation is just concatenating words, and reducing if 
necessary. 



Theorem V.2 (jOkada and Shudol (|200lh ) Let D 

and D' be two unfolded domains obtained by N times 
successive reflections of the same building block. If 
nD (l) — nD (7) f or an U sequence 7, then D and D' are 
transplantable, so also isospectral. 

Proof. Let G and G' — corresponding to D and £>', 
respectively — be the groups generated by the adjacency 
matrices: 



G=(tfW), G' = (N^); 



(101) 



then clearly G and G' are subgroups of the symmetric 
group Sat on N letters. Let F3 be the free group gener- 
ated by symbols o, b and c. Define the surjective homo- 
morphism 

$£> : F 3 H> G : 7 = 7172 . . . 7™ i-> M (7m) M (7m - l} . . . M (7l) . 

(102) 

Then 



G F 3 /ker$D and G' = F 3 /ker<f> D , 



(103) 



the latter notation being obvious. 

Now assume that n D {^f) — n D (7) for any sequence 7. 
Then 

ker$ D = {7 || ^(7) = 1} = {7 || n D ( 7 ) = N} 

= { 7 ||n 23 '( 7 )=^}={7l| ^'(7)=!} 

= ker$£,-. (104) 

(Note that $0(7) is a (0, l)-matrix, so that $^(7) = I if 
and only if n D (j) = N.) So the map 

A : Gh> G' : $13(7) i-> $73/(7) (105) 

yields an isomorphism between G and G'. 
Let (identity maps) 

p D : G^ GL(JV,C), p D ' : G' 1— > GL(iV, C) (106) 

be linear representations of G and G', respectively. Since 
the latter groups are isomorphic, 

p = p D ' o A : $ D ( 7 ) ^ $ D ,( 7 ) e GL(N, C) (107) 

is another linear representation of G. Since n ^) and 
n 15 (7) become (equal) characters of the representations 



2. Penrose-Lifshits mushrooms 

Since transplantation implies isolength spectrality, one 
might wonder if two billiards with the same length spec- 
trum are, in general, necessarily isospectral. 
M. Lifshits, exploitin g a constr uction attributed to 
R. Penrose (see, e.g., (|Rauchl . [l978l i). constructed a class 
of pairs of R 2 -domains that, while not isometric, have 
periodic geodesies of exactly the same lengths, including 
multiplicities. When the boundaries are (G°°) smooth, it 
follows that the two billiards have the same wave invari- 
ants, in the sense that the traces of their wav e groups , 
cosCt -y/A) , differ at most by a smooth function ( Melrose . 
19961 ). Such billiards provide drums that sound different 
but are similar geometrically. 

In this section we describe a construction of smooth 
Penrose-L ifshits mushroom pairs that ar e not isospectral, 



following ([Fulling and Kuchmenti . 120051 ) . The domains 
are smooth, so the spectral difference is not attributable 
to diffraction from corners. 

We start from a half-ellipse E with foci F and F' as shown 
in Fig. [15j The map 




F O F' 



FIG. 15 Starting half-ellipse. 



(109) 



whether applied to regions, curves, or points, indicates 
reflection through the minor axis of the ellipse. If objects 
are interchanged by that reflection, we call them dual. 
Now replace a line segment by a bounded smooth curve 
defined over the same interval, B\, on the left and B2 on 
the right, with B[ ^ B2 , to form a smooth domain f2 
(Fig. IT5|) . Finally, carry out the same replacement oper- 
ation (not self-dually) between the foci in two dual ways 
(M and M') to get two domains fii and H2 (Figs. ITTl 
and[j"8|). We call domains fii and constructed in this 
manner Penrose-Lifshits mushroom pairs, respectively. 



Theorem V.3 (jFuIling and Kuchmenti (|2005h ) // 

B\ and B2 are given and not dual, then there exist dual 
bumps M and M such that the resulting Penrose-Lifshits 
mushrooms Clj have the same length spectra and wave 
invariants but are not isospectral. 
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FIG. 16 Half-ellipse with two bumps. 



FIG. 18 Same as Fig. 1171 but perturbed in a dual way. 




FIG. 17 Perturbed half-ellipse with two bumps. 



Proof. Following ( Fulling and Kuchmentl . 2005 ) we first 
handle the length spectra ( MelroseT 1996 ; Zelditch . 
l2004bh . Geodesies i n an ellipse (logically) fall into 
two disjoint categor ies (|Berrvl . ll981tlKeller and Rubinowl . 
119601 : lRauchlll978h : those that intersect the major axis 



between the foci, and those that do so at or beyond the 
foci. (The smoothness assumption guarantees that the 
major axis will not bifurcate in flj by diffraction.) 
A similar observation holds for the domains Qj just de- 
scribed before Theorem IV. 31 Any geodesic originating in 
a curve B\ or B^ can never reach a curve M or M', and 
vice versa. 

The geodesies that do not intersect the focal segment 
FF' are the same for the two domains. Those for fii 
that do intersect this segment are identified one-to-one 
with their duals in ^2 by the reflection. Hence the two 
billiards are length isospectral. It now suffices to show 
that it is possible to choose an M in such a way that the 
two billiards are nonisospectral. One considers the lowest 
eigenvalue of the domains whose boundary is modified 
by a small perturbation. Using the Rayleigh-Hadamard 
formula for change i n the spectrum unde r domain 



perturbation s (cf. (Garabcdian and Schiffcr, 1952 



Ivanov et all Il977t IZelditchl l2004bl) or (jGarabedianl 



1998L section 15.1, Exercise 9)), one can prove that this 



eigenvalue is different for the two domains for a given 
choice of the perturbation. Thus, one can construct non- 
isospectral billiards having the same length spectrum, 
and the theorem is proved. ■ 



J. Analytic domains 

As all known counterexamples to the question "Can one 
hear the shape of a drum?" are plane domains with cor- 
ners, it might be possible that analytic drumheads are 



spectrally determine d. The paper of Zelditch ( 20091 ) is 



2004allcT ))" devoted to the 



part of a series (cf. (jZelditch 
inverse spectral problem for simply connected analytic 
Euclidean plane domains f2, the motivating problem be- 
ing whether generic analytic Euclidean drumh eads are 
determ ined by their spectra. The main results of lZelditch 
( 2009h give the strongest evidence to date for this conjec- 
ture by proving it for two classes of analytic drumheads: 
those with an up/down symmetry, and those with dihe- 
dral symmetry. 

Planar drumheads with symmetry. We now state 
the results more precisely. As before, by Lsp(f2) we de- 
note the length spectrum of f2, that is, the set of lengths 
of closed trajectories of its billiard flow. A bouncing ball 
orbit 7 is a two-leg periodic trajectory that intersects 
<9£1 orthogonally at both boundary points. By rotating 
and translating Q we may assume that 7 is vertical, with 
endpoints at A = (0,L/2) and B = (0, -L/2). Zelditch's 
iitchl . " 



inverse results ( (jZelditc 



2009)) pertain to the follow- 
ing two classes of drumheads, 2?i,l and X> m ,Z) which are 
defined as follows: (i) the class of drumheads with 
one symmetry a and a bouncing ball orbit of length 2L 
which is reversed by cr, and (ii) the class T> m ,L for m > 2 
of drumheads admitting the dihedral group D2 m (acting 
on m letters) as symmetry group and an invariant m-leg 
reflecting ray. 

Note that the class consists of simply connected 

real-analytic plane domains f2 with the property that 
there is an isometric involution a of f2 which "reverses" a 
nondegenerate bouncing ball orbit (that is, (7(7) = 7~ , 
i.e. the same orbit reversed) of lengt h = 2L. Ot her 
geometric properties can be found in (|Zelditchl . [2009l ). 



Let Spec B denote the spectrum of the Laplacian Ab 
of the domain Q with boundary conditions B. The 
result of Zelditch is that for Dirichlet (or Neumann) 
boundary conditions B, the map Spec B : 2?i,l h-> is 
one-to-one. As a corolla ry, o ne obtains the main r esult of 
llantchenko et all ([20021 ) and IZelditchl (|l999l [2000) that a 
simply connected analytic domain with the symmetries 
of an ellipse and with one axis of a prescribed length L 
is spectrally determined within this class. The above 
theorem admits a generalization to the special piecewise 
analytic mirror symmetric domains with corners that are 
form ed by reflecting the graph of an analytic function, 
see (Zelditch, 



„ I2009D . For m > 2, th at is, "dihedrally 
symmetric domains" , Zelditch! " (|2009l ) similarly proves 



that the map Spec s : r D m ,L 



5iV 



is one-to-one. 
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Higher dimensional drumheads with the 
symmetry of an ellipsoid. More generally, 



Hezari and Zelditch (2009) proved that bounded analytic 
domains Q C K n with ± mirror symmetries across all 
coordinate axes and with one axis height fixed (and 
also satisfying some generic non-degeneracy conditions) 
are spectrally determined among other such domains. 
That is, you can hear the shape of a real analytic drum 
in any number of dimensions if you know in advance 
that the mystery drums have the symmetries of an 
ellipsoid. It is one of the first positive higher dimensional 
inverse spectral results for Euclidean domains that is not 
restricted to balls. 



VI. EXPERIMENTAL AND NUMERICAL 
INVESTIGATIONS 

Although isospectrality is proved on mathematical 
grounds, the knowledge of exact eigenvalues and eigen- 
functions can not be obtained analytically for such sys- 
tems. Experimental as well as numerical simulations oc- 
curred very early in the history of billiards. In 1909, 
in the Bulletin international de l'Academie des Sci- 
ences de Cracovie, Zaremba proposed a way of "numer- 
ically" calculating solutions of th e Dirichle t and Neu- 

19091) 



mann problem at a given point (jZaremba . 1909). To 
solve the eigenvalue problem for the Helmholtz equation, 
one standard method is the bounda r y elem ent method 
(|Berrv and Wilkinsonl Il984t iRiddell . Il979h . However, 
this approach faces problems when the billiard h as cor- 
ners. Such situations have been addressed e.g. by IPisani 
(|l996h and lOkada et all (|2005bl ). The usual numerical 
methods to compute eigenvalues and eigenfunctions in 
polygonal billiards are based on the so-called "method 
of particular solu tions" introduce d by Fox, Henrici and 
Moler (FHM) in (|Fox et all 1 19671) . At a diffracting cor- 
ner with angle wa, a wave function \& admits a "corner" 
decomposition into Bessel functions valid at a distance 
smaller than to the nearest diffracting corner. In polar 
coordinates centered around the corner ira this decom- 
position reads 



tf(r, 6) 



k 



cikJka(kr) sin(fca6>), 



(110) 



where J„ are Bessel functions of the first kind, and 
k = \f~E. The sine function in Eq. (|110j) ensures that 
the function ^(r, 9) is zero on the boundary edges con- 
nected to corner ira. The idea of FHM is to require that 
^ also vanish on the rest of the boundary at a finite num- 
ber of points, and to truncate the sum (|110[) . This gives 
a system of m linear equations, which admits a non-zero 
solution {ctfc, 1 < k < m}if and only if the matrix cor- 
responding to this linear system is singular. The FHM 
method therefore consists in varying the energy E and 
tracking the singularities of the matrix M. 
Unfortunately, for more than one diffracting corner it be- 
comes virtually impossible to track singularities, espe- 
cially since in various circumstances the FHM method 
fails to converge when the number of terms included in 
Eq. (|110j) is increased. Even for the paradigmatic pair 
with half-square base shape (Fig. |TJ), which is one of the 



simplest isospectral billiards, each pair has four diffrac- 
tive angles: two 37r/2 and two 37r/4 angles, and the FHM 
method fails to give eigenvalues with a good accuracy. 
This is why attention has been focused on physical ex- 
periments. 

All known pairs of isospectral billiards are built on the 
same principle as the "historical" pair 7% of Fig. LTJ As 
explained in section |H] any initial building block pos- 
sessing three sides along which to unfold the block can 
be used to construct an isospectral pair. In particular, 
the properties of the resulting pair will depend strongly 
on the choice of the initial building block. Physicists 
have mainly concentrated on the paradigmatic example 
of Fig. [T] This allows us to make comparisons between 
the different approaches. 

In this section we review both experimental and numeri- 
cal investigations which give insight into the behavior of 
eigenvalues and eigenfunctions for isospectral billiards. 



A. Numerical investigations 

1. Mode-matching method 

Numerical approaches to the study of isospectrality for 
the billiards of Fig. [TJ have followed the experiments of 
Sridhar and Kudrolli that will be reviewed in the next 
subsection. Various approaches have been used in or- 
der to solve the Helmholtz equation (A + E)^ = with 
Dirichlet boundary conditions inside the billiards. The 
first numerical results were obtained by Wu , Sprung and 
Martorell and reported in (|Wu et all Il995l ). Using the 



mode- matching method described in section IV.F. 21 they 
obtained eigenvalues of the billiard as the values for which 
the determinant of the matrix M, given by Eq. (|95j) . van- 
ishes. The results obtained by this method are displayed 
in column 2 of Table [TTl As expected, both billiards yield 
the same values. The numerica l resul ts were found to 
vary linearly in 1/N. IWu et all (|l995h compared their 
results to results obtained by a finite-difference method 
consisting in discretizing the Laplacian A. This finite- 
difference method gives the results displayed in column 
1 of Table [TT] (the numerical results are again exactly the 
same for both billiards) . As a check of the validity of this 
approach, one can identify the eigenvalues of triangular 
states. The lowest-energy triangular state are expected 
to have eigenenergies equal to 57r 2 /d 2 and 107r 2 /d 2 . As 
one can see in Table [n] these eigenvalues correspond to 
the ninth and 21st mode, respectively, consistently with 
Table U 



2. Expansion of eigenfunctions around the corners with the 
domain-decomposition method 

The main draw back of the mode-matching method of 
(|Wu et all Il995l) is the fact that one has to know analytic 
solution s of the Helmho ltz equation on subdomains of the 
billiard. Driscolll ( 1997) used a numerica l meth od based 
on an algorithm by Descloux and Tollevl ( 19831 ). partic- 
ularly suited to treating the case of polygonal billiards. 
The idea is again to decompose the billiard into domains, 
each domain T>i containing only one diffracting angle a%. 
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On each domain the restriction of the eigenfunction ^/ 
is supposed to be some \fj that admits a Bessel func- 
tion expansion around corner a%, according to Eq. (|110[) . 
Truncation of this expansion to some finite order reduces 
the problem to that of finding the coefficients of the ex- 
pansion for the ^i. Mode-matchin g numerically leads 
to und esired singularities. Instead, iDescloux and Tollev 
(|1983T ) used an algorithm minimizing a function that mea- 
sures discrepancies between the ^ and between their 
derivatives at the boundaries between subdomains. Im- 
provement of this algorithm allowed Driscoll to obtain the 
first 25 eigenvalues for both billiards of Fig, [pwith an ac - 
curacy of up to 12 digits. iBetcke and Trefethenl (|2005h 
used a modified method of particular solutions using 140 
expansion terms at each singular corner, 140 boundary 
points on each side of the polygon, and 50 interior points 
to obtain following estimates for the first three eigenval- 
ues: 2.537943999798, 3.65550971352 and 5.17555935622. 



B. Experimental realizations 

1. Electromagnetic waves in metallic cavities 

Many experimental studies have been carried out 
on chaotic quantum billiards to check the various 
properties conj e cture d analytically for chaotic systems 



(jBohigas et all . 119841 ). One commonly used method 



is based on the correspondence between the stationary 
Schrodinger equation and the Hclmholtz equation for 
electromagnetic waves in two dimensions (which is also 
the equation obeyed by vibrating plates). The exper- 
iments are carried out by sending electromagnetic mi- 
crowaves into a cylindrical copper cavity. The height 
h of the cavity is small, and the two other dimensions 
are shaped according to the desired billiards to investi- 
gate. For wavelengths A > 2h, i.e. frequencies below 
vq = c/2h, all modes obey the two-dimensional wave 
equation (A + k 2 )^> = 0. The E z component of the elec- 
tric field plays the role of the quantum wave and vanishes 
on the boundary. Probes allow one to send an electromag- 
netic wave into the cavity and to measure the transmis- 
sion spectrum. In particular, eigenvalues correspond to 
resonances in the transmission spectrum. Various choices 
of the probe locations ensure that no resonance is missed. 
If A is the area of the cavity, the number of resonances 
below Vq is approximately given by -^(ir/h) 2 . But the 
quality factor of the cavity is proportional to h\ therefore 
one has to find a compromise between a high quality fac- 
tor and a large number of resonances. 
Measurements of the intensity of the wave function (or 
here the electri c field) were a c hieve d by the perturbation 
body method (jSridhar et all Il992h : The resonance fre- 
quency of the cavity is shifted by the presence of a small 
metallic body inside the cavity. This shift is a function of 
the square of the electric field at the point of the metallic 
body. 

The first experimental investigation of isospectral bil- 
liar ds was realized at N o rtheas tern University, Boston, 
by ISridhar and Kudrolli (1994). Sridhar and cowork- 
ers carried out various studies on chaotic quantum 
billiards, such as the Sinai billiard (a square billiard 
with a circular obstacle in the interior), and the Buni- 



movitch stadium- shaped billiard , observing the scarring 
of ei genfunction s (Sridhar, 1991) or localization phenom- 
ena (|Sridhar and Helleit ll992) for such billiards. The ex- 
periments aimed at investigating isospectrality were real- 
ized on cavities having the shape of the isospectral pair 
of Fig. [[J 

Experimentally, each cavity has nine rectangular sides. 
The base shape is an isosceles rectangular triangle (a half- 
square) whose smaller side is d — 76 mm (3 in.) long. The 
height of the cavity is h = 6.3 mm (~ 0.25 in.), so that 
microwaves at frequencies below vq = 25 GHz are actu- 
ally two-dimensional. Measurements carried out to ob- 
tain the 54 lowest eigenvalues showed that, as expected, 
the eigenvalues of the two cavities are equal. Relative 
discrepancies of 0.01 to 0.2% between pairs of eigenval- 
ues were found. These discrepancies and the width of 
the resonances were assumed to be caused by imperfec- 
tions due to the assembly of the pieces forming the cavity. 
This experiment also allowed insight into the properties 
of eigenvalues of isospectral pairs. It was checked that 
the eigenvalues found experimentally agree with the Weyl 
formula (|5T|) for the integrated density of states: 

77(E) ~ £-E- ^-VE + IC. (Ill) 

For the choice d = 3 in. one gets an area A = 31.5 in. 2 
and a perimeter C = 27 in.; the constant JC is given by 
Eq. (j53"j) and yields JC = 5/12. It was observed that at 
least for the lowest eigenvalues no degeneracy occurred. 
By measuring the electric field inside the cavity, some of 
the lowest eigenfunctions were obtained. The results for 
the ten first eigenvalues are displayed in Table HU It is 
interesting to note that these pairs of eigenfunctions look 
quite different, although they possess the same eigenval- 
ues. It was checked that one eigenfunction could be de- 
duced from the other by transplantation. The particular 
case of the ninth mode, which is a triangular state, is 
well reproduced. Indeed, as shown in Table [TTl the mea- 
sured 9-th eigenvalue is very close to its theoretical value 

E = b TT 2 /d 2 . 

Later [Dhar et aL ( 2003 ) applied a similar technique to 
a chaotic isospectral billiard made of the billiard with 
half-square base tile with scattering circular disks inside, 
showing experimentally that isospectrality is indeed re- 
tained, provided scatterers are added in a way consistent 
with the unfolding rules. 



2. Transverse vibrations in vacuum for liquid crystal smectic 
films 



Another experimental realization of Kac's membranes 
was achieved using liquid crystal films in a smectic phase, 
spann ing a shape of the form of the isospectral bil- 
liard ([Even and Pieranski . Il999h . First, the shapes were 
etched in circular stainless-steel wafers of diameter 4cm 
and thickness 125 /im. The smectic film is then drawn on 
the shape, and after a few hours it reaches an equilibrium 
with uniform thickness e of several hundred nanometers 
(corresponding to a few dozens monomolecular layers) 
over the whole surface. The whole experiment is set in 
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vacuum. The film then obeys the wave equation 

d 2 z 



jAz = pe 



dt 2 ' 



(112) 



where 7 is the intrinsic tension of the film (in the ex- 
periments 7 ~ 5.10~ 2 N/m), and p is the density, with a 
vertical displacement z vanishing on the border. The film 
is excited by a voltage applied by an electrode under the 
film, and the amplitude and phase of its oscillations are 
measured by sending in a laser beam and measuring its 
deviations with a photodiode. The signal detected is pro- 
portional to the height of the film at the position of the 
electrode. The frequency of the excitation is varied from 
a few Hertz to several kiloHertz, and eigenfrequencies cor- 
respond to resonance peaks. Displacing the electrode over 
the whole shape allows to reconstruct eigenmodes. 
The experiment was carried out on isospectral billiards 
with an isosceles triangular base shape: two angles /3 = 7 
are equal, while the third one is varied from a = 67.5° 
to 97.5°. The angle a — 90° corresponds to the exam- 
ple of Fig. [TJ The first 30 modes for both shapes were 
measured. The average relative difference between two 
eigenvalues for a given mode is 0.3%, which is within the 
estimated experimental error of order 0.5%. For the right 
angle triangle (a = 90°) the modes can be compared with 
other numerical or experimental results. Data for the ten 
first eigenvalues of t he a = 90° billiards were given by 
Even and Pier anski (1999) and are displayed in Table El 



When the parameter a is varied, there is an avoided cross- 
ing between eigenvalues of the eighth and ninth mode. 
Since the ninth mode is a triangular mode (see section 
IV.F.1[) and the eighth is not, the coupling between these 
two modes necessarily comes from experimental imper- 
fections. 

This experiment has also been tested on a billiard where 
the gluing scheme of the base triangles is modified. That 
is, tile E in Fig. [13] (left) is flipped around the line 
x = 5d/2. This leads to a significantly different spec- 
trum. In particular, the triangular modes are no longer 
eigenstates of such a billiard. Again this is a check that 
the way the tiles are glued together, according to the rules 
constructed from finite projective spaces or from Sunada 
triples, is of primary importance for isospectrality. 



3. Isospectral electronic nanostructures 



Recently (jMoon et aD,[2008) an experiment was done in- 
volving electrons confined in isospectral billiards, with 
the purpose of using transplantation to reconstruct the 
quantum phase of measured wavefunctions. Each bil- 
liard consisted of a wall of 90 CO molecules, constructed 
by positioning the molecules with the tip of a scanning 
tunneling microscope. The chosen billiards were built ac- 
cording to the pattern of Fig. [TJ but the base shape was 
cho sen to be a triangle with angles (tt/2,%/3, 7r/6). As 
in ( Even and Pieranski . 19991) . it was checked that bil- 
liards violating the isospectral construction rule led to a 
different res ult. 

Amusingly, M oon et all (120081 ) took Kac's question lit- 
erally by converting the average measured spectra into 



audio frequencies, checking that one could indeed "hear" 
non- isospectrality. 



VII. SUNADA THEORY 

The examples of isospectral billiards considered so far 
can be proved to be isospectral by quite simple tools. 
However historically they were constructed by a group- 
theoretical approach. The mathematical theory of 
isospectrality rests on a theory by Sunada. We first re- 
view the necessary basic notions of group theory. Then, 
in section lYlI.El we introduce Sunada Theory. 



A. Permutations 

Following the usual conventions, we denote permutation 
action exponentially (i.e. the image of an element x by 
the permutation g is x 9 ) and let elements act on the right. 
We denote the identity element of a group by id or 1, if no 
special symbol has been introduced for it before. A group 
G without its identity id is denoted G x . The number of 
elements of a group G is denoted by \G\. A permutation 
group (G, X) is a pair consisting of a group G and a set 
X such that each element g of G defines a permutation 
g : X — > X of X, and the permutation defined by the 
product gh, g,h <E G, is given by gh : X — > X : x 1— > 
(x 9 ) h . An involution in a group is an element g of order 
2, that is, such that g 2 = id. 



B. Commutator notions 

The group-theoretic setting of Sunada theory requires 
introduction of some notions such as the commutator of 
two groups and perfect groups. The conjugate of g by h 
is g h = h~ 1 gh. Let H be a group. The commutator of 
two group elements g, h is equal to [g, h] = g~ 1 h~ 1 gh. 
The commutator of two subsets A and B of a group G 
is the subgroup [A, B] generated by all elements [a, 6], 
with a £ A and i e B. The commutator subgroup of G 
is [G, G], also denoted by G' . Two subgroups A and B 
centralize each other if [A, B] = {id}. The subgroup A 
normalizes B if B a = B for all a € A, which is equivalent 
to [A, B] being a subgroup of B. HA and B are two 
subgroups of the group G, then they are conjugate(d) 
if there is an element g of G such that A 9 — B. The 
subgroup A of G is (a) normal (subgroup) in (of) G if 
A 9 = A for all g £ G. In such a case, we write A<G. If 
A ^ G, we also write A < G. 

Inductively, we define the nth central derivative [G, G]r„i 
of a group G as [G, [G, G] an d the nth normal 
derivative [G,G] (n) as [[G, G\( n -i)t [ G > G ](n-i)]- For 
n = 0, the zeroth central and normal derivatives are 
by definition equal to G itself. If, for some natural 
number n, [G,G] (n) = {id}, and [G,G] ( „_ 1} ^ {id}, 
then we say that G is solvable (soluble) of length n. If 
[G, G][ n ] = {id} and [G, G][ n _x] ^ {id}, then we say that 
G is nilpotent of class n. The center of a group is the 
set of elements that commute with every other element, 
i.e., Z{G) = {z G G || [z,g] = id,V 5 G G}. Clearly, if a 
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rank 


Finite differences 


Mode matching 


Electromaj 


;netic waves 


Smcctic fiims (reiative values) 


1 


1.028936 


1.028535 


1.02471 


1.02481 


1.000000 


1.000000 


2 


1.481865 


1.481467 


1.46899 


1.47194 


1.438000 


1.430000 


3 


2.098249 


2.097467 


2.08738 


2.08831 


2.040000 


2.027000 


4 


2.649715 


2.649547 


2.64079 


2.63985 


2.571000 


2.548000 


5 


2.938176 


2.937434 


2.93297 


2.92949 


2.854000 


2.823000 


6 


3.732689 


3.732334 


3.72695 


3.71892 


3.623000 


3.570000 


7 


4.295193 


4.294728 


4.28393 


4.28388 


4.184000 


4.153000 


8 


4.677665 


4.677532 


4.67021 


4.66917 


4.554000 


4.507000 


9 


5.000002 


5.000000 


4.98838 


4.98531 


4.861000 


4.811000 


10 


5.291475 


5.290275 


5.27908 


5.27278 


5.150000 


5.095000 


11 


5.801531 


5.801138 


5.78755 


5.78371 






12 


6.433894 


6.432156 


6.41357 


6.43781 






13 


6.866260 


6.866226 


6.84891 


6.84718 






14 


7.159802 


7.159343 


7.15242 


7.16045 






15 


7.694737 


7.692417 


7.67783 


7.70604 






16 


8.463655 


8.463257 


8.44285 


8.45947 






17 


8.613536 


8.611169 


8.57859 


8.62220 






18 


9.012405 


9.010349 


8.99495 


8.97209 






19 


9.609968 


9.609791 


9.60312 


9.59562 






20 


9.921131 


9.921040 


9.92583 


9.93689 






21 


10.000008 


10.000000 


10.00330 


10.03932 






22 


10.571020 


10.569736 


10.55227 


10.55740 






OQ 


1 i nftftm ft 


1 i nft^vov 
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24 


11.419551 


11.418850 


11.41874 


11.40569 






25 


11.984650 


11.984080 


11.99364 


11.98033 







TABLE II Comparison between the first eigenvalues Ei of the isospectral pair obtained by various methods, expressed in 
units of 7r 2 /d 2 . The ninth mode corresponds to the triangular mode: its normalized e igenvalue i s exp ected to be equal to 
5. (The conversion from frequencies to lengths is done assuming vacuum in the cavity; IWu et all (|l995l ) gave the values for 
electromagnetic cavities with a factor of 1.0006 corresponding to the presence of air in the cavity.) 



group G is nilpotent of class n, then the (n — l)th central 
derivative is a nontrivial subgroup of Z(G). 

A group G is the central product of its subgroups A and 
B if AB = G, A n B is contained in the center of G, 
and A and B centralize each other. Sometimes we write 
G = A o B in such a case. 

A group G is called perfect if G = [G, G] = G'. 

Let R be a finite group. The Frattini group 4>(R) of R is 
the intersection of all proper maximal subgroups, or is R 
if R has no such subgroups. 



C. Finite simple groups 

A group is simple if it does not contain nontrivial normal 
subgroups. 

The finite simple groups are often regarded as the elemen- 
tary particles in finite group theory. Before we explain 
this more precisely, recall that a composition series of a 
group G is a normal series 



1 = H < H x < ■ ■ ■ < H n = G, 



(113) 



such that each ifj is a maximal normal subgroup of Hi+\. 
Equivalently, a composition series is a normal series such 
that each factor group Hi+i/H is simple. The factor 
groups are called composition factors. 
A normal series is a composition series if and only if it 
is of maximal length. That is, there are no additional 
subgroups that can be "inserted" into a composition se- 
ries. The length n of the series is called the composition 
length. 

If a composition series exists for a group G, then any 
normal series of G can be refined to a composition series. 



Furthermore, every finite group has a composition series. 
A group may have more than one composition series. 
However, the Jordan-Holder theorem states that any two 
composition series of a given group are equivalent. 



The classification of finite simple groups (see ([Solomon . 
l200ll ) for a survey) states that every finite simple group is 
cyclic, or alternating, or is contained in one of 16 families 
of groups of Lie type (including the Tits group, which 
strictly speaking is not of Lie type), or one of 26 sporadic 
groups. 



Conwav et al. ( 19851) provided a li st of the finite simple 



groups, see also (|Gorensteinl . Il980l p. 490-491). In this 
review, we encounter several aspects of certain simple 
groups in the construction theory of counter examples to 
Kac's initial question. 



D. p-Groups and extra-special groups 

The present section will be useful for construction of ex- 
amples in section IVII.FI 

For a prime number p, a p-group is a group of order p n 
for some natural number n 7^ 0. A Sylow p-subgroup of 
a finite group G is a p-subgroup of order p n such that 
p n+1 does not divide \G\. 

A p-group P is special if either [P, P] — Z(P) — <fi(P) is 
elementary Abelian or P itself is. (A group is elementary 
Abelian if it is Abelian, and if there exists a prime p such 
that each of its nonidentity elements has order p.) Note 
that P/\P,P] is elementary Abelian in that case. So 



P/[P,P]^V(n,p), 



(114) 



where V(n,p) is the n-dimensional vector space over ¥ p 
(here seen as its additive group), and \P\ = p n \[P, P]\. 
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Hence we have the exact sequence 

1^[F,F]^P4 V(n,p) i-» 1. 



(115) 



If furthermore \Z(P)\ = \[P,P}\ = \<f>(P)\ = p, P is called 
extra-special. 

We now present a classification for extra-special groups 
that depends on the knowledge of the nonabelian p- 
groups of order p 3 . 

There are four no nabelian p-groups of order p 3 — see 
(|Gorensteinl . ll980l) . First we have M — M(p): 



M(p) = (x, y, z\\x p = y p = z p = 1, 

[x,z] = [y,z] = 1, [x,y] = z). (116) 

(Note that this is the general Heisenberg group of order 
p 3 which we will encounter later on.) Next, define 

M 3 (p) = (x,y\\ x? 2 =yP= l,x* - x^ 1 ). (117) 

Finally, we have the dihedral group D of order 8 and the 
generalized quaternion group Q of order 8. 



Theorem VII. 1 (Go renstein An extra- 

special p-group P is the central product of r > 1 
nonabelian subgroups of order p 3 . Moreover, we have the 
following. 



(1) If p is odd, P is isomorphic to N k M r k , while if 
p = 2, P is isomorphic to D k Q r ~ k for some k. In 
P 2 



either case, \P\ — ■" 2r + 1 



(2) If p is odd and k > 1, N k M r ~ k is isomorphic to 
NM 7 " -1 , the groups M r and NM r ~ l are not iso- 
morphic and M r is of exponent p. 

(3) If p — 2, then D k Q r ~ k is isomorphic to DQ r ~ l if 
k is odd and to Q r if k is even, and the groups Q r 
and DQ r ~ x are not isomorphic. 

(All the products considered are central products.) 



E. Sunada Theory 

We now turn to the main theorems of Komatsu and 
Sunada, which allowed Gordon et al. to produce the 
first known example of isospectral billiards. Sunada's 
idea was to reduce the problem of finding isospectral 
manifolds to a group-theoretical problem, namely, 
constructing triplets of groups having a certain property. 
As the groups that appear in Sunada's proof are Galois 
groups, we need some more definitions. 

A field extension L/K is called algebraic if every element 
of L is algebraic over K, i.e., if every element of L is a root 
of some non-zero polynomial with coefficients in K. (Field 
extensions which are not algebraic, i.e. which contain 
transcendental elements, are called transcendental.) 
Let K be an algebraic number field of degree n. Recall 
that a number field is a finite, algebraic field extension of 
Q; its degree is the dimension over Q as a Q- vector space. 
A standard example is 



The ring of integers of an algebraic number field K, often 
denoted by Ok, is the ring of algebraic integers contained 
in K. An algebraic integer is an element of K that is a 
root of some monic polynomial with coefficients in Z. 
The (Dedekind) zeta function Qk{s) (associated with K), 
s being a complex variable, is defined by 



(118) 



taken over all ideals I of the ring of integers Ok of K, 
/ ^ {0}. Note that N%(I) denotes the norm of I (to Q), 
equal to \Ok/I\- 

An ideal P of a ring R is a prime ideal if it is a proper 
ideal and if for any two ideals A and B in R such that 
AB C P, we have that A C P or B C P. Let p be a 
rational prime. Let P\ , . . . , P g be the prime ideals of Ok 
lying above p. Then 



where 



e l = e K (Pi)- 



(119) 



(120) 



Here e^iPi) is the ramification index of Pi over K. If 
ej > 1 for some i € {1, . . . , g}, then p is said to be 
ramified in K. If ej = 1 for all i, p is unramified in K. 

The conjugate elements of an algebraic element a, over 
a field K, are the roots of the minimal polynomial of a 
over K. (For example, the cubic roots of 1 are 1, — 1/2 + 
\/3/2i, —1/2 — v3/2i. The latter two roots are conjugate 
elements in the field K = Q[\/— 3].) 

Let K = Q(9) be as above, that is, an algebraic number 
field of degree n (9 G C). Suppose 9i, 62, ■ ■ ■ , n are the 
conjugates of 9 over Q. If 



7i) = ■ ■ ■ = \ayu n ) = K, 
then K is a Galois extension of Q. 



(121) 



Suppose that E is an extension of the field F (written 
as E/F). Consider the set of all automorphisms of E/F 
(that is, isomorphisms a from E to itself such that 
a(x) = x for every x € F). This set of automorphisms 
with the operation of function composition forms a 
group, sometimes denoted by Aut(E/F). If E/F is a 
Galois extension, then Aut(E/F) is called the Galois 
group of (the extension) E over F, and is usually denoted 
by Gal(E/F). 

A number-theoretic exercise which asks for non- 
isomorphic number fields Ki and K2 with the same zeta 
function has the following answer: 



Theorem VII.2 (jKomatsul fll976D) Let K be a finite 
Galois extension o/Q with Galois group G = Gal(K/Q), 
and let Ki and K2 be the subfields of K corresponding 
to subgroups G\ and G2 of G, respectively. Then the 
following conditions are equivalent: 

(i) Each conjugacy class of G meets G\ and G2 in the 
same number of elements; 
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(ii) The same primes p are ramified in Ki and K2 and 
for the unramified p the decomposition of p in ¥L\ 
and K2 is the same; 

(iii) The zeta functions of Ki and K2 are the same. 

In particular, if G\ and G2 are not conjugate in G, then 
and K 2 are not isomorphic while having the same 
zeta function. It should be noted that several such triples 
(G, Gi , G2) are known — see the examples in this section. 
Any group triple (G, Gi, G2) satisfying Theorem lVII.2f i) 
is said to satisfy "Property (*)" . 

Sunada's idea was to establish a counterpart of this the- 
orem for Riemannian geometry. In that context, there 
is an analogue for the Dedekind zeta function. For CDt a 
Riemannian manifold, one defines 



Cot(s)=X>7 s , 9te(s)»0, 



where 



< Ai < A 2 < 



(122) 



(123) 



are the non-zero eigenvalues of the Laplacian for 971. The 
function (jm has an analytic continuation to the whole 
plane, and it is well-known that C<xfi 1 (s) = £ajt 2 (s) if and 
only if SUli and 9Jt2 are isospectral. 

The following theorem gives sufficient conditions for two 
manifolds to have the same zeta function. 



Theorem VII.3 (jSunadal <|l985h ) Let n : DJl i-> 9Jt 

be a normal finite Riemannian covering with covering 
transformation group G, and let tt± : 9Jli i-> 9Jto and 
7T2 : WI2 1 y 9Jto be the coverings corresponding to the 
subgroups Hi and H2 of G, respectively. If the triplet 
(G, -Hi, #2) satisfies Property (*), then the zeta functions 
Csuti ( s ) an d Can 2 ( s ) are "identical. 

The proof of the latter theorem makes use of an interest- 
ing trace formula, which we present now. 
If A is a non-negative self-adjoint operator of a Hilbert 
space, one defines the trace of A as an extended real num- 
ber by the possibly divergent sum ^2 k {Aek,ek), where 
{ej}j is an orthonormal base of the space. It is of trace 
class if and only if Tr(A) < 00. 

Let V be a Hilbert space on which a finite group G acts 
as unitary transformations and let A : V H> V be a 
self-adjoint operator of trace class such that A commutes 
with the G-action. For a subgroup H of G, denote by 
V H the subspace of if-invariant vectors. 

Trace Formula. The restriction of A to the subspace 
V G is also of trace class, and 



tv(A\ v a 1 ) = tx(A\ v a 2 ). 



(125) 



Even if G\ and G2 are not conjugate, the manifolds 9Jti 
and 9Jl 2 could possibly be isometric. 



Theorem VII. 4 (jSunadal (|l985T There exist finite 
coverings 7i"i : 9Jti M> 9Jto and ir 2 : M2 n- 9Jto of Riemann 
surfaces with genus > 2 such that for a generic metric 
go on 9Ho, the surfaces (9Jti,7r 1 g , o) an d (Ct^^So) are 
isospectral, but not isometric. 



Sunada's theorem allows us to construct isospectral pairs 
provided we find triples (G, Gi,G2) satisfying Property 
(*) — "Sunada triples". 
Now we give examples of such triples. 



F. Examples of Sunada triples 



Example 1 — see iGerstl (|l970h . Let G be the 

semidirect product Z/8Z X k Z/8Z, and define Gi and 
G 2 by 



d ={(1,0), (3,0), (5,0), (7,0)}, 
G 2 ={(1,0), (3,4), (5,4), (7,0)}. 



(126) 



Example 2 — see iGassmann (192 



Let G 



be the symmetric group on six letters {a, b, c, d, e, /}. Set 
Gi = {1, (06) (cd), (ac)(bd), (ad){bc)} (127) 

and 

G 2 = {1, (ab)(cd), (ab)(ef), (cd)(ef)}. (128) 
Example 3 — see iKomatsul ( 197fjh . Let G 2 and G 2 



be two finite groups with the same order, and suppose 
that their exponents (equal to the least common multiples 
of the orders of their elements) both equal the same odd 
prime p. Set |Gi| = |G 2 | = p h for h € N x and embed 
Gi and G2 in the symmetric group S p h on p h letters by 
their left action on themselves. For a conjugacy class [g] 
corresponding to the partition 



P + jjH +p, 



we have 



\([g}nG 1 )\=p h -l = \([g]nG 2 )\, 



(129) 



(130) 



while \([g] (1 G,)| = otherwise. 

Concretely, let Gi = (Z/pZ) 3 , and let G2 be the group 



tv(A\ v a) = Yl (\G g \)-hv(gA), (124) 
[ff]e[G] 

where [G] = {[3]}, [g] is the conjugacy class of g in G 
and G g is the centralizer of g in G. 

If the triplet (G, Gi,G2) satisfies Property (*), then 



G 2 = (a, Ml a p = b p = [a,b] p = 1, 

a[a, b] = [a,b]a,b[a,b] = [a, 6] 6), (131) 

that is, G2 is the extra-special group of order p 3 . Then 
(S p 3,Gi,G 2 ) verifies Property (*). 

One can in fact generalize Komatsu's example by defining 
the following group. The general Heisenberg group H„ of 
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dimension 2n + 1 over F g , with n a natural number, is 
the group of square (n + 2) x (n + 2)-matrices with entries 
in ¥q, of the following form (and with the usual matrix 
multiplication) : 



(132) 



where a, (3 £ F™, c € F 9 , and with 1„ the n x n-unit 
matrix. Let a, a' , /3, /3' G F™ and c, c' € ¥ q ; then 






Q 






a' 




• 


In 






I„ 



















1 a + a' c + c' + (a,/3') 
= 10 In P + 

1 



(133) 



Here (x,y), with a; = (x\,X2, ■ ■ ■ ,x n ) and y = 
(2/i,j/2, • ■ • , 2/«) elements of F™, denotes £12/1+222/2 + • ■ • + 

Xnyn • 

The following properties hold for H„. 

(i) H n has exponent p if q = p with p an odd prime; 
it has exponent 4 if q is even. 

(ii) The center of H„ is given by 



been concentrated on mixed Dirichlet-Neumann bound- 
ary conditions, that is, having either = or d n tp = 
on different intervals of the boundary (n being the nor- 
mal to the boundary). This is much simpler than the 
Dirichlet problem. Simple instances of m ixed-boundary 



cond ition isospe ctral pairs are proposed in (jLevitin et al. 
20061 ) (see also ( Jakobson et al. . 20061 )). Their simplest 
example is reproduced in Fig. Q3B The eigenfunctions are 
given by 

7r(m + l/2)x imy 
sin-^ 1 ' sin— p n>l,m>0, (135) 



for the square of size d, and 

tt(to + 1/2):e . n(n + l/2)j/ 
sin = sin = 

7r(n + l/2)x 7t(to + 1/2)2/ 
— sin t= sin = , 

for the triangle of size dy/2. 



(136) 

m > n > 0, 




{(0,c,0) || ceFj. 
(iii) H„ is nilpotent of class 2. 



(134) 



Then, as above, (S p 2n+i , H„, (Z/pZ) 2,l+1 ) verifies Prop- 
erty (*). 

Any finite group arises as the fundamental group of a 
compact smooth manifold of dimension 4. For a triplet 
(G, Gi, G2) of the type described in Example 3, we find 
a compact manifold DJIq with fundamental group G. Let 
071 be the universal covering of 2Jl - Then the quotients 
DJli = VJl/Gi have non- isomorphic fundamental groups 
G h i = 1,2. By Theorem EECl the manifolds (Wti, it*g Q ) 
and (9JI2, Tr^ffo) are isospectral for any metric go on SEJto, 
but not isometric. 



VIII. RELATED QUESTIONS 

The literature on isospectrality is large, and it is out of 
the question to review the entire field. In the present pa- 
per we have concentrated on the questions addressed by 
planar two-dimensional domains with Dirichlet boundary 
conditions. To open the topic further, we now mention 
some questions related to the main one discussed in the 
present paper, some of which have been addressed in the 
literature, and some of which remain open problems. 



A. Boundary conditions 

So far we have mainly dealt with billiards with Dirich- 
let boundary conditions. More recently attention has 



FIG. 19 Isospectral billiards with mixed Neumann-Dirichlet 
boundary conditions. Solid line, Dirichlet; dashed line, Neu- 
mann. 



These examples can be generalized: Levitin et al. ( 20061) 
gave a procedure to construct similar pairs. The idea is to 
construct an elementary domain, or "construction block" , 
whose boundary is made of two line segments a and b 
on the plane, with ends joined by two arbitrary curves. 
Imposing any mixed Neumann-Dirichlet boundary con- 
ditions on the construction block boundary, one obtains 
a Neumann-Dirichlet isospectral pair by gluing the con- 
struction block together with its reflection with respect 
to either a or b (and imposing Neumann boundary condi- 
tions to the segment itself, Dirichlet to its image). This 
technique can be further generalized by gluing together 
more copies of the construction block, yielding more com- 
plicated examples. In particular this method shows that 
for mixed boundary conditions it is possible to construct 
isospectral pairs such that one member is connected and 
the other is not; isospectral pairs such that one member 
is smooth and the other is not; isospectral 4-tuples; bil- 
liards whose spectrum remains invariant when Dirichlet 
and Neumann bound aries are swapp ed. These billiards 



were investigated by iJakobson et 



(120061) . The sim- 
plest example is a billiard of semi-circular shape: if the 
equation of the billiard on the complex plane is given by 
{z G C;0 < arg(z) < n; \r\ < 1}, the Dirichlet boundary 
conditions correspond to {z £ C; |r| = 1, 7r/4 < arg(z) < 
37r/4, !ft(z) < 0}. A necessary condition for this Dirichlet- 
Neumann isospectrality is that the Dirichlet boundary 
has the same total length as the Neumann boundary. 
Such domains have been investigated numerically 
([Driscoll and Gottlieb! . 120031 ) as well as analytically 
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(jOkada and Shudol. 120011). and experimental s etups have 
been proposed bv lDriscoll and Gottlieb! (|2003l ). 
All these examples have the property that the length dif- 
ference between the Dirichlet boundary and the Neumann 
boundary are the same. This turns out to be a necessary 
condition similar to those obtained from Wevl's law ( | 52 | ) 



applyi ng to isospectral billiards derived by iLevitin et al 



(|2006l ) for mixed-boundary condition isospectral billiards. 
In particular, such isospectral pairs need to have the 
same area, the same length difference between the Dirich- 
let boundary and the Neumann boundary, and the same 
curvature-singularity properties, namely, the quantity 



2 / n{s)ds+y2 ~ 

JdB DD 







NN 



DN 



2f3 2 



(137) 

where n is the curvature and j3 represents the angles at 
the Dirichlet-Dirichlet, Neumann-Neumann or Dirichlet- 
Neumann boundary intersections, must be the same for 
both billiards. 

Finally, w e observe tha t, for s ome of the examples pro- 
duced by ILevitin et al. ( 2006t ), it was shown that two 
isospectral domains produce a different number of nodal 
domains (domains separated by nodal lines where \f r = 0; 
see section IV.H|) . 



B. Homophonic pairs 

Homophonic pairs in K 2 are nonisometric compact 
domains that have a distinguished point such that the 
corresponding (normalized) Dirichlet cigenfunctions take 
equal values at that point. This could be interpreted in 
the following way: If the corresponding drums are struck 
at these special points, then they sound the same in 
such a way that every frequency is excited to the same 
intensity for each. 

An example of two billiards t hat are isospectral and 
homophonic ( Buser et all , 1994 ) is provided in Appendix 
lAl (example 21 1 right). These billiards sound the same 
when struck at the interior points where six triangles 
meet. 



C. Spectral problems for Lie geometries 

There exists a vast literature on spectral prob- 
lems for (finite) graph s - - see the excellent paper 
( van Dam and Haemersl 120031 ). In this section we con- 
sider a spectral ("Kac type") problem for graphs that are 
associated with the most important incidence geometries. 
We have seen in the previous sections that the construc- 
tion of isospectral pairs is based on properties of finite 
projective spaces and their automorphism groups. In this 
section we show that this construction is a special case of 
a wider class of similar constructions based on so-called 
generalized polygons, which are the natural generaliza- 
tion of projective planes. 

One defines a finite axiomatic projective plane II of order 
n, where n € N, as a point-line incidence structure satis- 
fying the following conditions: (i) each point is incident 
with n + 1 lines and each line is incident with n + 1 points; 



(ii) any two distinct lines intersect in exactly one point 
and any two distinct points lie on exactly one line. One 
also traditionally requires that n be > 2 to exclude the 
uninteresting cases of a single line and a point not on it 
(n = — 1), a single line and one point on it (n = 0), or the 
three vertices and three sides of a triangle (n = 1). This 
is equivalent to requiring that II contains an ordinary 
quadrangle (four points with no three on a line) as sub- 
geometry. It is easily seen that a finite projective plane 



of order n has 



1 points and n 



1 lines. 



The obvious examples of finite projective planes are 
the projective planes PG(2, q) over finite fields ¥ q as 



defined in section |IVl In this case the order n = \¥ q \ is a 
prime power, and in fact no examples of finite projective 
planes of non prime power order are known. A classical 
theorem of Moufang states that a finite projective plane 
is isomorphic to some PG(2,g) if and only if a certain 
configurational property corresponding to the classical 
theorem of Desargues is satisfied. Projective planes of 
this type are therefore often called Desarguesian, and 
since these correspond to planes coordinatized over finite 
fields, we also use this terminology for projective spaces 
of dimension n > 3, as already mentioned. However, 
many finite projec tive planes are known w hich are not 
Desarguesian; see (jHughes and Piperl . Il973h . 



Generalized polygons. Let n > 3 be a natural number. 
A (thick) generalized n-gon or (thick) generalized polygon 
(GP) is a point-line geometry T = (V,B,T), where V is 
the point set, B is the line set and I C (V x B) U (B x V) 
is a symmetric incidence relation, so that the following 
axioms are satisfied: 

(i) r contains no fc-gon (in the ordinary sense) for 2 < 
k < n; 

(ii) Any two elements x, y G V U B are contained in 
some ordinary n-gon in T; 

(iii) There exists an ordinary (n + l)-gon in T. 

The point graph of a point-line geometry is the graph of 
which the vertices are the points of the geometry, and 
for which two vertices are joined by an edge if they are 
collincar in the geometry. Equivalently, a generalized 
polygon could be defined as a point-line geometry for 
which the point graph is bipartite of diameter n and girth 
2n (see, e.g., Fig. EH}. 

The generalized 3-gons are precisely the aforementioned 
projective planes. If (iii) is not satisfied for T, then 
r is called thin. Otherwise, it is called thick. Each 
thick generalized n-gon, n > 3, T has an order, there 
are (not necessarily finite) constants s > 1 and t > 1 
so that each point is incident with t + 1 lines and 
each line is incident with s + 1 points. We then say 
that r has order (s,i). Note that, for a point x and a 
line L, xIL means that (x, L) € I (and so also (L, x) € I). 

Collinearity matrices and a spectral problem. 

Suppose r = (VjB,!) is a finite GP (r has a finite 
number of points and lines) of order (s,i), set IT^I = v, 
and let {x\, X2, ■ • ■ , x v } = V be the point set. Define the 
collinearity matrix A(r) — A as the v x w-matrix (Oy) 



for which a, 



1 if Xi 



7^ xi (the latter notation 
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FIG. 20 The unique generalized quadrangle of order 2: the 
symplectic W(3, 2). 



meaning that Xi and different collinear points), 

and otherwise. So it is the adjacency matrix of the 
point graph of T. The (point) spectrum of T is the 
spectrum of A, and we denote it by spec(A). 

The following quantum mechanical question is the Kac 
inverse problem for the theory of GPs. 

Question VIII. 1 Let T and V be distinct finite thick 
generalized polygons with associated collinearity matrices 
A and A', respectively. Does spec(A) = spec(A') imply 

that r = r ? 

Clearly, a similar problem can be posed for the line spec- 
trum, but as points and lines play essentially the same 
role in a GP, we only consider the question in its above 
form. 

Question IVlTLTI can be reduced to an important question 
in the theory of GPs: 



Theorem VIII.2 (jThasI (|2007ah ) Let r and V be dis- 
tinct finite thick generalized polygons with associated 
collinearity matrices A and A! , respectively. Then 
spec (A) = spec (A') if and only if V and V have the 
same order. 

Details of the proof can be found in Appendix |B] 



D. Further questions 

As mentioned, the literature on isospectrality is huge and 
continuously growing. There is also a vast literature on 
isospectral graphs. In this section we state some funda- 
mental open problems (which reflect our personal inter- 
est) on billiards and graphs. 

Interesting problems in construction theory are numer- 
ous: We state only some of them. Perhaps the single 
most important open problem in Kac theory is the fol- 
lowing: We have constructed pairs of isospectral billiards 
made of 7, 13, 15, or 21 tiles. Is it possible to go beyond 
that number? In mathematical words, can one show that 
for all N € N there exists an N* > N such that there are 
isospectral pairs on N* tiles? Equivalently, can one show 
that there are infinitely many pairs of involution graphs 
that yield isospectral pairs? 



All examples constructed so far are polygonal examples. 
Even if different base tiles can be chosen, the unfold- 
ing rule imposes the presence of corners in the bound- 
ary of the billiard. A natural question is thus whether 
one can construct isospectral ]R 2 -domains with smooth 
boundaries. 

We have seen that point-line duality in finite projective 
spaces is at the root of billiard isospectrality and provides 
a transplantation property between billiards. Since only 
one recipe is known for constructing isospectral pairs, one 
may ask the following: Is it possible to construct isospec- 
tral pairs which are not transplantable? More generally, 
are the following statements achievable: Derive criterions 
for pairs of involution graphs to yield isospectral plane do- 
mains; Construct isospectral pairs on oo tiles (perhaps by 
a free construction); Find examples of (planar) isospec- 
tral pairs not coming from Sunada triples, or still arising 
from Sunada triples but not being transplantable. 

On the group theoretical level, we pose the following ques- 
tion: Are the operator groups of (transplantable) isospec- 
tral pairs always two-transitive? If so, the classification of 
finite simple groups could be used to classify such opera- 
tor groups. In the same spirit, one could ask as to whether 
other finite simple groups can act as operator groups. A 
related question is to develop a theory of isospectral "do- 
mains" on general buildings. Note that the proj ective 
completion of Mr is a rank 2-building over R (see (|Titd . 
1974 ) for an introduction on buildings). The same ques- 
tions could all be formulated for "isospectral n-tuples", 
n > 2. 



ACKNOWLEDGMENTS 

K. T. acknowledges the Fund for Scientific Research - 
Flanders (Belgium) for financial support. This paper 
was partly written while K. T. was visiting the Discrete 
Mathematics group of the Technical University of Eind- 
hoven, The Netherlands, whose hospitality he gratefully 
acknowledges. The paper was finished while both au- 
thors were hosted by Institut Henri Poincare (IHP) at 
Paris, whose hospitality is gratefully acknowledged. 



Appendix A: Gallery of examples 
1. Some modes 

Here we plot some eigenfunctions for the pair of billiards 
of Fig. [TJ Figure [2T] corresponds to the fundamental 
mode, Fig. [22] to the first triangular mode, whose nodal 
lines coincide with edges between the triangular tiles. 
Figure l23l corresponds to an excited state. 



2. The 17 families of isospectral pairs and their 
mathematical construction 

The following gallery presents the 17 k nown fami- 
lies o f isospect r al pair s, as obtained by iBuser et ~al\ 
(Il994h : iGiraudl lj2005bh : and lOkada and Shudol (l200lh . 
All are based on a Sunada triple (G.Gi.G^), where 
G = PSL(n + l,q) is the special linear automorphism 
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FIG. 21 (Color online) Fundamental mode. 




FIG. 22 (Color online) First triangular mode (ninth mode). 



group of a finite projective space of (q n+1 — l)/(q — 1) 
points, and Gi,G2 are two subgroups, generated by 
ax,bi,c\ and a2,&2,C2 given below, respectively. These 
automorphisms are collineations of order 2 of the un- 
derlying finite projective space; ai,i>i, ci act on points 
while a2,&2,C2 act on hyperplanes, numbered from to 
(q n+1 — l)/(q — 1) — 1. The generators a, and &i allow to 
construct the graphs (see section III.Aj) that specify the 
way in which the tiles are glued together. Figures [Ml to 
[40] give examples of pairs of isospectral billiards obtained 
by applying the unfolding rules on an equilateral triangle 
(left panel) or on a scalene triangle (right panel). 

Interestingly, the structure of pairs 136 and 152 forbids 
the construction of any proper billiard, that is, struc- 
tures where triangles do not overlap. It is quite simple 
to convince oneself of this fact. In the case of the bil- 
liard 136 (see Fig. [32]) . the initial triangle is unfolded six 
times around each of its corner. Clearly, to have a non- 
overlapping billiard, each angle should be less than 7r/3, 
which is impossible unless the initial triangle is equilat- 
eral. 

For the billiard 152 (see Fig. [37]) . the initial triangle is 
unfolded six times around two of its corners, four times 
around the third one, and thus two angles have to be less 
than 7r/3 and one less than tt/2. While it is possible to 
construct such a billiard, it is impossible to get a pair 
of planar billiards. Indeed, the role of the angles is ex- 
changed from one billiard to the other, which leads to the 
condition that the three angles be less than tt/3. On the 
other hand, the presence of a loop in the pair 21 1 requires 
that one angle of the base triangle be tt/3. 



FIG. 23 (Color online) 15th mode. 




FIG. 24 Pair 7i. Sunada triple G = PSL(3,2), Gi = 

{ai,bi,a), i = 1,2, with ai = (0 1)(2 5), &i = (0 2) (3 4), ci = 
(0 4)(1 6), a 2 = (0 4) (2 3), b 2 = (0 1)(4 6), c 2 = (0 2)(1 5). 




FIG. 25 Pair 7 2 . Sunada triple G = PSL(3, 2), Gi = 

{<H,bi,Ci), i = 1.2, with oi = (0 1)(2 5), 6i = (1 5) (3 4), ci = 
(0 4)(1 6), a 2 = (0 4)(2 3), b 2 = (0 6)(1 4), c 2 = (0 2)(1 5). 




FIG. 26 Pair 7 3 . Sunada triple G = PSL(3, 2), Gi = 

(oj,6i,Ci>, i = 1,2, with oi = (2 5) (4 6), 6i = (1 5) (3 4), ci = 
(0 4)(1 6), a 2 = (0 3) (2 4), b 2 = (0 6)(1 4), c 2 = (0 2)(1 5). 




FIG. 27 Pair 13i. Sunada triple G = PSL(3,3), d = 

(<H,bi,Ci), i = 1,2, with oi = (0 12)(1 10) (3 5) (6 7), 
bi = (0 10)(2 9) (3 4)(5 8), ci = (0 4)(1 6) (2 11)(9 12), 
a 2 = (0 4) (2 3) (6 8) (9 10), b 2 = (0 12)(1 4) (5 11) (6 9), 
c 2 = (0 10)(1 5)(2 7) (3 12). 
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FIG. 28 Pair 13 2 . Sunada triple G = PSL(3, 3), d = 



(<H, bi,Ci), i = 1,2, with ai = 
bi = (1 12)(2 9)(3 8) (4 5), ci 
a 2 = (0 4)(2 3) (6 8) (9 10), b 2 
c 2 = (0 10)(1 5)(2 7)(3 12). 



(0 12) (1 10) (3 5) (6 7), 
(0 4)(1 6) (2 11) (9 12), 
(0 1)(4 12)(5 11)(8 10), 






FIG. 30 Pair 13 4 . Sunada triple G = PSL(3,3), d = 
{ai,bi,a), i = 1,2, with oi = (1 7)(3 5) (4 9) (6 10), 
6i = (0 5)(1 2)(6 12)(9 11), ci = (0 4)(1 6)(2 11) (9 12), 
a 2 = (0 9) (4 10)(6 8)(7 12), b 2 = (0 11)(1 8) (2 7)(3 4), 
c 2 = (0 10)(1 5)(2 7)(3 12). 






FIG. 29 Pair 13 3 . Sunada triple G = PSL(3,3), Gi = 

(ai,bi,Ci), i = 1,2, with oi = (1 7)(3 5) (4 9) (6 10), 
bi = (1 12) (2 9) (3 8) (4 5), ci = (0 4)(1 6) (2 11)(9 12), 
a 2 = (0 9)(4 10) (6 8)(7 12), b 2 = (0 1)(4 12) (5 11) (8 10), 
c 2 = (0 10)(1 5) (2 7) (3 12). 



Appendix B: Spectral problems for Lie geometries 



1. Generalized polygons 



Generalized polygons were introduced by iTits (1959) 
in order to have a geometric interpretation of certain 
Chevalley groups of rank 2. They are also the building 
bricks of (Tits) buildings, the natural geometries for the 
groups with a BN-pair. 

A group G is said to have a BN-pair (B,N), where 
B, N are subgroups of G, if the following properties are 
satisfied: (BN1) (B,N) = G; (BN2) H = B n N < N 
and N/H = W is a Coxeter group (see, e.g., (|Titsl . 
19741) ) with distinct generators s±, S2, • ■ • , s n ; (BN3) 
Bs.iBwB C BwB U BsiwB whenever to € W 
and i € {1, 2, . . . , ra}; (BN4) Sl B Si ^ B for all 
i € {1,2,..., n}. The subgroup B, respectively W, is a 
Borel subgroup, respectively the Weyl group, of G. The 
natural number n is called the rank of the BN-pair. 

Example. Suppose PG(1, q) is the projective line over 
the finite field F^; so PG(l,g) has q+l points. Consider 
the natural action of PSL(2,g) on PG(l,q), and let 
x and y be distinct points of the projective line. Set 
B = PSL(2,g) 2 , and N = PSL(2, q){ x , v y Then (B,N) 
is a BN-pair for PSL(2, q). Here N/(b'(~) N) = W is just 
the group of order 2. 

Example. Consider the Desarguesian projective plane 
PG(2,q), and PSL(3, q) in its natural action on the 
latter plane. Let (x, L) be an incident point-line pair, 
and A a triangle (in the ordinary sense) that contains x 
as a point and L as a side. Set B = PSL(3, q)r x ,L) an< i 
N = PSL(3, q) a ; then (B, N) is a BN-pair for PSL(3, q) 
and N/ (B n N) — W is the dihedral group of order 6. 



See dPavne and Thasl. Il984t iThas et all 120061: iThasi 







FIG. 31 Pair 13 5 . Sunada triple G = PSL(3,3), d = 
(ai,bi,Ci}, i = 1,2, with oi = (1 7)(3 5)(4 9)(6 10), 
bi = (0 5)(1 2) (6 12) (9 11), ci = (0 4)(1 6)(2 11) (9 12), 
a 2 = (0 9)(4 10) (6 8) (7 12), b 2 = (0 11)(1 8)(2 7)(3 4), 
c 2 = (0 10)(1 5)(2 7) (3 12). 



consider only thick CPs. 

Now let G be a group with a BN-pair (B,N) of rank 2. 
One can associate a generalized polygon B(G) with the 
group G in the following way. For this purpose, define 
Pi = (B,B S1 ) and P 2 = (B,B S2 ). 

• Call the right cosets of Pi "points" . 

• Call the right cosets of Pi "lines" . 

• Call two such (distinct) cosets "incident" if their 
intersection is nonempty (so P\ g is incident with 
P 2 h, g,heG, if P ig C\P 2 h^%). 

Then B(G) is a CP — say a generalized n-gon for some 
natural number n — on which G acts naturally as an au- 
tomorphism group that permutes transitively the ordered 
n-gons (in the ordinary sense) . 

Conjecture [J. Tits, (|TitsL Il974 
§11.5.1)]. If a finite thick generalized 
n-gon is such that the automorphism group 
permutes transitively the ordered n-gons (that 
is, if A is associated with a BN-pair), then 
A is isomorphic with the GP of an absolutely 
simple group over a finite field, or with the 



GP of a Ree group of type 2 F4 over a finite 
field. 



For more on the clas sificati on of BN-pairs of rank 2 . 
see ([Thasl . 120091 12010[) and (|Thas and Van Maldeghem 
2008). 



2004 Van Maldeghem . 1998 ) for standard references on 
the subject of generalized polygons. In this paper we 



2. Duality principle 

Let T = (V,B,I) be a GP of order (s,t). Then T D = 
(B,P,1) clearly again is a GP, but now of order (t,s). 
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FIG. 32 Pair 13 6 . Sunada triple G = PSL(3, 3), d = 

{ai,h,Ci), i = 1,2, with oi = (0 2)(1 7) (3 6) (5 10), 
6i = (0 6)(2 4) (3 8) (5 9), a = (0 5)(1 2) (6 12) (9 11), 
a 2 = (0 7) (3 11)(6 8)(9 12), b 2 = (0 8)(1 10) (5 11)(7 9), 
c 2 = (0 11)(1 8) (2 7)(3 4). 




FIG. 34 Pair 13 8 . Sunada triple G = PSL(3,3), d = 
(ai,bi,Ci), i = 1,2, with oi = (0 10) (1 5)(2 7)(3 12), 
bi = (0 4)(2 3)(6 8) (9 10), ci = (0 5)(1 2) (6 12) (9 11), 
a 2 = (0 4)(1 6) (2 11)(9 12), b 2 = (0 12) (1 10) (3 5) (6 7), 
c 2 = (0 11)(1 8)(2 7) (3 4). 



\2 



7VV 



7 




FIG. 33 Pair 13 7 . Sunada triple G = PSL(3, 3), d = 

{ai,bi,a}, i = 1,2, with oi = (0 2)(1 7)(3 6) (5 10), 
h = (0 4) (2 3) (6 8) (9 10), ci = (0 5)(1 2) (6 12) (9 11), 
a 2 = (0 7) (3 11) (6 8) (9 12), b 2 = (0 12) (1 10) (3 5) (6 7), 
c 2 = (0 11)(1 8) (2 7) (3 4). 



(The latter geometry is called the point-line dual of L.) 
So any theorem which holds for a GP, has a dual inter- 
pretation; we call this principle "duality principle" . 



3. Automorphisms and isomorphisms 

Let T = (P,B,I) and V = (V',B', V) be GPs. Then an 
isomorphism between T and V is a pair (a, j3) for which a 
is a bijection between V and V' , /? is a bijection between 
B and B', and zLL (in T) if and only if x a lL^ . If there 
is an isomorphism between T and T', we say they are 
"isomorphic", and write r = V. 

If r = r' one speaks of an "automorphism". The set 
of all automorphisms of a GP forms a group, and the 
classical examples of GPs are those examples that are 
associated with a Chevalley group (or, equival e ntly, with 
a "classical" BN-pair), cf. ([Van Maldegheml . Il998ri for 
more details. 



4. Point spectra and order 

Let r be a finite thick GP of order (s, t), with associated 
collinearity matrix A. Our first concern is to calculate 
spec(A). 



First, we recall the theorem of iFeit and Higmanl (|1964T ): 
A finite thick generalized n-gon exists if and only if 
n G {3,4,6,8}. We will do a case-by-case analysis 
according to this result. 

Case n = 3. 

Recall that a generalized 3-gon is the same as an ax- 
iomatic projective plane. Now let T be a finite projective 
plane of order n, n > 2, and put n 2 + n + 1 = 0, its num- 
ber of points. Then A = J„ — I D , where J„ is the all 1 
x 0-matrix, and I„ the x 0-identity matrix. It follows 
that 

spec(A) = {-1,0 - 1}. (Bl) 
So if the spectra of two finite projective planes coincide, 





FIG. 35 Pair 13 9 . Sunada triple G = PSL(3,3), d = 
(oi,6i,c<), i = 1,2, with oi = (0 10) (1 5) (2 7) (3 12), 
bi = (1 10) (3 6)(5 7) (9 11), ci = (0 5)(1 2)(6 12) (9 11), 
a 2 = (0 4)(1 6)(2 11)(9 12), b 2 = (0 3) (2 4)(6 8)(7 11), 
c 2 = (0 11)(1 8)(2 7) (3 4). 



their orders do as well. 
Case n = 4. 

Let L be a thick generalized 4-gon, or also "generalized 
q uadrangle" (GQ) o f orde r (s, t). Then using the results 
of IPavne and Thasl (|l984 1.2.2), we have 



spec(A) = {-t - 1, s - 1, s(t + 1)}. 



(B2) 



Now let r' also be a thick GQ, with the same spectrum, 
of order (s',t'). There is only one negative eigenvalue, 
so -t - 1 = -t' - 1 and t = if. Since s - 1 < s(t + 1) 
(s' — 1 < s'(t' + 1)), it also follows that s = s', and hence 
r and r' have the same order. 

Case n = 6. 

For this case, we need one more definition. A distance 
regular graph G with diameter d is a regular connected 
graph with valency k for which there exist natural num- 
bers bo = k, bi, . . . , bd-i', ci = 1, C2, . . . , e<2 such that for 
each pair of vertices x and y at distance j, we have 

. \r j _ 1 (y)nr 1 (x)\=c j ,l<3 <d; 

. \r j+1 ( y )nr 1 (x)\ = b j ,o<j<d-i. 

Now define the co nstants a, = |r,fa) D Fi(a;)| for < 



j < d. A result of iBrouwer et all (| 19891 ) claims that the 
eigenvalues of the point graph of G are the eigenvalues of 
the following " intersection matrix" : 



/0 1 

k a\ C2 
bi a 2 



\ 



V 



Crf 

bd-i a d J 



(B3) 



It just so happens to be that the collinearity graph of 
a generalized hexagon (a generalized 6-gon) is distance 
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FIG. 36 Pair 15i. Sunada triple G = 

PSL(4,2), d = (ai,bi,*), i = 1,2, with 
ai = (0 14) (1 12) (2 6) (4 5)(7 11)(9 10), 6i = 
(1 13)(2 7)(4 6)(8 9), ci = (1 14)(2 12)(3 4)(8 11), a 2 = 
(0 11)(1 5)(3 4)(6 10)(8 9)(13 14), b 2 = (0 10)(1 2)(6 9)(12 14), 
c 2 = (0 5) (2 4) (6 7) (11 14). 






FIG. 37 Pair 15 2 . Sunada triple G = 

PSL(4,2), Gi = ((H,bi,Ci), i = 1,2, with 
at = (0 14)(1 12) (2 6) (4 5) (7 11)(9 10), bi = 
(1 13)(2 7) (4 6) (8 9), ci = (0 12) (1 6) (3 5) (7 8), a 2 = 
(0 11)(1 5) (3 4)(6 10) (8 9)(13 14), b 2 = (0 10) (1 2) (6 9) (12 14), 
c 2 = (0 13)(1 11)(2 3) (7 10). 



FIG. 38 
PSL(4,2), 



with 



Pair 153. Sunada triple G 

G t = (ai,bi,Ci), i = 1,2, 
ai = (0 14) (2 11)(4 7)(5 6) (8 10) (12 13), 6i = 
(1 13) (2 7)(4 6) (8 9), ci = (0 12) (1 6) (3 5) (7 8), a 2 = 
(0 9)(2 5) (3 4) (6 8) (10 11)(12 13), b 2 = (0 10) (1 2)(6 9)(12 14), 
c 2 = (0 13) (1 11)(2 3) (7 10). 





FIG. 39 Pair 15 4 . Sunada triple G = 

PSL(4,2), d = (a^b^a), i = 1,2, with 
oi = (0 14) (2 11)(4 7)(5 6) (8 10) (12 13), 6i = 
(1 4)(2 8)(7 9)(8 9), ci = (0 12) (1 6)(3 5)(7 8), a 2 = 
(0 9)(2 5)(3 4) (6 8)(10 11)(12 13), b 2 = (6 9) (7 13)(12 14), 
c 2 = (0 13)(1 11)(2 3) (7 10). 



regular with a diameter of 3. An easy exercise yields, for 
a thick generalized hexagon of order (s,t), the following 
intersection matrix: 



B 



/ 1 \ 

s(t+ 1) s- 1 1 

st s - 1 f + 1 

\ (i + l)(s - 1) / 



(B4) 



The determinant of B — zLj has the following roots: 

x = — t— 1, at = s(t+l), se = s — 1 — Vsi, x = s— l+Vst. 

(B5) 

One observes that — t — 1 is strictly the smallest eigen- 
value, while s(t + 1) is the largest. It now easily follows 
that if a generalized hexagon of order (s',t') has the 
same spectrum as T, then it has the same order. 

Case n = 8. Let V be a thick generalized octagon (8-gon) 
of order (s, t). Again, the point graph is distance regular, 
now with a diameter of 4. The intersection matrix is 
easily seen to be the following: 



/ 





1 













s(t + 1) 


s- 1 


1 













st 


s - 1 


1 













st 


s - 1 


t + 1 


V 











st (t - 





B 



which has eigenvalues 

-t- 1, 

s - 1 - 



x 



x 



'1st, 



= s(t- 

l + « 



-1), 

/ 2si. 



\ 



(B6) 



(B7) 



The third largest eigenvalue is s — 1, so if V is a thick 
generalized octagon of order (s',t') with the same spec- 
trum, then s = s'. As s(t + 1) is the largest eigenvalue 



of spec (A), it follows that t = t'. This ends the proof of 
Theorem |WL2l ■ 



5. Concluding remarks 

In this section, we make some comments on generalized 
polygons that are characterized by their order. 

Projective planes. For some small values, e.g. 
n = 2, it is known that there is a unique projective 
plane of order n (up to isomorphism). It is well-known, 
however, that as soon as n is large enough and not a 
prime, nonisomorphic examples exist. On the other 
hand, for p a prime, only one example is known, namely 
the classical example PG(2,p) arising from a BN-pair in 
PSL(2,p). 

Generalized quadrangles. Many infinite classes of 
GQs are known, and several examples with small param- 
eters are completely det ermined by their order. We refer 
the interested reader to ( Payne and Thasl 1984 . Chapter 
6) for these examples. We make some comments accord- 
ing to the known orders. Below, q is always a prime 
power. We also assume s < t by re asons of du ality. [De- 
tails and references can be found in (|Thasl . l2004L Chapter 
3).] 

• (s, t) = (g 2 ,<7 3 ). Only one example is known 
(for each q), namely, the Hermitian quadrangle 
H(4,q 2 ). 

• (s,t) = (q — l,q + 1). If q > 8 and q is even, 
nonisomorphic examples are known for every q. In 
the other cases, only unique examples are known. 

• (s, t) = (q, q). If g is odd, nonisomorphic examples 
are known for every q. If q > 8 and q is even, we 



40 




FIG. 40 
PSL(3,4) 



triple G 
= 1,2, 



with 



ai = (2 7) (3 11)(5 12)(8 18)(13 14)(15 17) (16 20) 

fci = (0 17) (3 8) (4 12) (6 13) (9 19) (14 15) (16 18) 

ci = (1 8) (2 16) (4 11)(5 19) (7 14) (10 17)(13 20) 

a 2 = (0 1)(4 17) (7 12) (9 16)(10 20)(11 13) (15 19) 

b 2 = (0 20) (3 16) (6 11) (8 15) (9 19) (10 12) (14 18) 
c 2 = (1 8) (2 16) (4 11)(5 19) (7 14) (10 17) (13 20). 



have the same remark. The other values give unique 
examples. 

• (s,t) — (q,q 2 ). If q > 5, nonisomorphic examples 
are known for every q. The examples of order (2, 4) 
and (3, 9) are unique. 

Generalized hexagons. Up to duality, only two 
classes of generalized hexagons are known (both asso- 
ciated with classical groups): the split Cayley hexagons 
H(g>) of order q, q a prime power, and t he twisted trialit 



hexag ons T(q,q 3 ) of order (g,q ), ci. ([Van Maldeghe: 



•My 
3ml . 



119981 . Chapter 2). We k now that H(q) ^ H(q) D if and 
only if q is a power of 3 (|Van Maldeghem . 19981) . If q is 
not a power of 3, H(g) ^ H(q) D , while both have the 
same spectrum. 

Generalized octagons. Up to duality, the only 
known thick finite generalized octagons are the Ree-Tits 
octagons 0(q), where q is an odd power of 2; they can 
be co nstructed from a BN-pa ir in the Ree groups of type 
2 F 4 (|Van Maldegheml . Il998l Chapter 2). They have 
order (q, q 2 ). 



of lengths of closed geodesies 7j . 

Denote by (T*M, J2j dxj A d£j) the cotangent bundle of 
M equipped with its natural symplectic form. Given the 
metric g, we define the metric Hamiltonian by 



n(x,o = \z\ = 



n+l 



Y^^i-nUr (C2) 



and define the energy surface to be the unit sphere bundle 

S*M = {(x,Z)\\\Z\ = l}. (C3) 

The geodesic flow G* is the Hamiltonian flow 

G* = exptEff : T*M \0^ T*M \ Q, (C4) 

where 3# is the Hamiltonian vector field. Since it is 
homogeneous of degree 1 with respect to the dilatation 
(x,£) H> (x,r£), r > 0, one can restrict G l to S^M. Its 
generator is also denoted by S. 

Livsic's cohomological problem asks whether a cocycle 
/ G C°°(S*M) satisfying 



fds = 



(C5) 



for every closed geodesic of the metric g is necessarily a 
coboundary / = where 5 is the generator of the 

geodesic flow G and g is a function with a certain degree 
of regularity. Under a deformation g e of a metric g = go 
preserving the extended Lsp(M, g) (including multiplici- 
ties), one has 



L 



gds = 0, VA. 



(C6) 



When the cohomology is trivial, one can therefore write 
g = S(/) for some / with the given regularity. One 
does not expect the cohomology to be trivial in general 
settings, but the results might be interesting for the 
length spectral deformation problem. 



Final ly, see 
1998) for 



more 



the monograph 
information on 



( Cvetkovic et al 



graph spectra. 

van Dam and Haemerd |2003) surveyed the known 
cases of graphs that are determined by their spectrum. 
Some generalized quadrangles with small parameters 
are mentioned that are uniquely determined by their 
spectrum. Since such examples must h ave the property 
that t hey are determined by their order. iPavne and Thas 
( 19841 . Chapter 6) also yielded these examples. 



Appendix C: Livsic cohomology 

In this appendix we describe a connection between 
isospectrality and cohomology. Let (M, g) be a Rieman- 
nian manifold without boundaries. The length spectrum 
is the discrete set 



Lsp(M,g) = {L 



7i 



< L 72 < 



} 



(CI) 
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